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ABSTRACT 


Various  methods  of  constructing  a  set  of  mutually  orthogonal  latin  squares 
are  presented  and  the  theoretical  aspects  of  various  methods  are  discussed. 
Illustrative  examples  of  constructing  latin  squares  and  sets  of  mutually  ortho¬ 
gonal  latin  squares  are  given.  The  methods  of  constructing  latin  squares  and 
sets  of  orthogonal  latin  squares  are  complete  and  partial  confounding,  frac¬ 
tional  replication,  analysis  of  variance,  group,  projecting  diagonals,  ortho¬ 
morphism,  pairwise  balanced  design,  oval,  code,  product  composition,  and 
sum  composition.  The  methods  of  construction  designated  as  partial  confound¬ 
ing,  fractional  replication,  analysis  of  variance,  and  sum  composition  appear 
not  to  have  been  discussed  previously  in  the  literature.  The  methods  of 
complete  confounding  and  of  projecting  diagonals  have  been  discussed;  the 
actual  construction  procedure  has  been  illustrated  with  several  examples. 

The  sum  composition  method  has  interesting  consequences  in  combinatorial 
theory  as  well  as  in  the  construction  of  orthogonal  latin  squares.  Lastly, 
equivalences  of  fourteen  combinatorial  systems  to  orthogonality  in  latin  squares 
has  been  investigated  and  described. 


SOME  TECHNIQUES  FOR  CONSTRUCTING  MUTUALLY 
ORTHOGONAL  LATIN  SQUARES 

12  3 

W.  T.  Federer  ,  A.  Hedayat  ,  E.  T.  Parker 

4  5  6 

B.  L.  Raktoe  ,  Esther  Seiden  ,  and  R.  J.  Turyn 


I.  Introduction  and  Some  Terminology 

The  purpose  of  this  paper  Is  to  present  a  set  of  methods  for  constructing 
mutually  orthogonal  latin  squa  es  and  to  exhibit  some  squares  produced  by  each 
of  the  methods.  The  set  of  methods  presented  herein  was  discussed  in  a  series 
of  informal  seminars  held  during  the  weeks  of  July  14-18  and  21-25,  1969,  by  the 
authors  at  Cornell  University.  The  motivation  for  these  discussion  was  derived 
from  results  obtained  by  Hedayat  [1969]  and  from  the  optimism  of  the  authors. 

New  procedures  for  constructing  a  set  of  mutually  orthogonal  latin  squares  and 
new  views  of  present  methods  of  construction  were  desired  in  order  to  advance  the 
theory  of  mutual  orthogonality  in  latin  squares. 


Professor  of  Biological  Statistics,  Cornell  University  and  Visiting  Professor, 
Mathematics  Research  Center,  University  of  Wisconsin  (on  sabbatical  leave 
1969-70). 

2 

Assistant  Professor,  Cornell  University. 

3 

Professor  of  Mathematics,  University  of  Illinois,  and  Visiting  Professor, 
Cornell  University  (July,  1969). 

4 

Associate  Professor,  University  of  Guelph  and  Visiting  Associate  Professor, 
Cornell  University  (January  to  August,  1969). 

5 

Professor,  Michigan  State  University,  and  Visiting  Professor,  Cornell  Uni¬ 
versity  (June,  July,  August,  I969). 

6 

Mathematician,  Raytheon  Corporation,  and  Visiting  Professor,  Cornell  Uni¬ 
versity  (July,  1969). 


Sponsored  by  the  United  States  Army  under  Contract  No. :  DA-31-124 -ARO-D- 
462. 


As  may  be  noted  from  the  table  of  contents,  the  different  sections  were 
written  by  different  authors.  An  attempt  was  made  to  have  a  consistent  notation 
and  a  uniform  style.  Although  much  more  work  is  required  to  finalize  the  method 
in  several  of  the  sections  enough  is  known  about  the  method  to  use  it  to  con¬ 
struct  a  latin  square  of  any  order  or  to  construct  a  set  of  two  or  more  mutually 
orthogonal  latin  squares.  Also,  a  number  of  equivalences  may  be  noted  for  some 
of  the  methods. 

The  theory  of  mutual  orthogonality  in  latin  squares  has  application  in  the 
construction  of  many  classes  of  experiment  designs  and  in  many  combinatorial 
systems.  The  latter  subject  is  discussed  in  section  XV  where  the  equivalences 
of  various  combinatorial  systems  are  presented.  With  regard  to  the  former  sub¬ 
ject,  there  is  an  ever  present  need  for  new  experiment  designs  for  new  experi¬ 
mental  situations  in  order  for  the  experimenter  not  to  have  to  conduct  his  ex¬ 
periment  to  fit  known  experiment  designs. 

Some  of  the  notation  and  terminology  that  will  be  utilized  is  presented 

below. 

Definition  1. 1 .  A  latin  square  of  order  n  on  a  set  with  n  distinct  elements 
is  an  n  X  n  matrix  each  of  whose  rows  and  columns  is  a  permutation  of  the  set 

l  ■ 

Example: 

is  a  latin  square  of  order  3  on  =  {1,2,3}  • 
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Definition  L  2  .  Two  la  tin  squares  L.  =  (a,.)  and  L_  =  (b . . )  of  order  n  are 
-  1  lj  2  ij 

,  ,  2 

said  to  be  orthogonal  if  the  n  ordered  p.lirs  (a,, ,  b, , )  (i,j  =  l,2,...,n)  are 
all  distinct.  Note  that  and  need  not  be  defined  on  the  same  set. 
Example; 

and 

Definition  I.  3.  The  members  of  a  set  of  t  latin  squares  Lp  L  ,  . ..  ,  L  of 
order  n  are  said  to  be  mutually  (pairwise)  orthogonal  if  L.  is  orthogonal  to 
Lj,  i#j,  i,  j  =  1,2,. . .  ,t  .  Hereafter  by  an  0(n,t)  set  we  mean  a  set  con¬ 
sisting  of  t  mutually  orthogonal  latin  squares  of  on.  T  n  . 

Example; 

,  and 


Latin  squares  and  orthogonal  latin  squares  have  at  least  187  years  of  history. 
Hedayat  [1969],  Section  IX  has  presented  a  reasonably  good  picture  of  this 
history  which  will  not  be  repeated  here.  It  is  planned  to  prepare  a  historical 
account  of  developments  related  to  orthogonality  in  latin  squares  and  to  publish 
this  material  together  with  a  bibliography  elsewhere. 
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II.  Factorial  Confounding  Construction  of  O(n.t)  Sets 
II.  1.  Complete  Confounding 

A  factorial  treatment  design  consists  of  all  possible  combinations  of  two 

i 

or  more  factors  each  at  two  or  more  levels.  The  set  of  all  combinations  of  m 
factors  each  at  n  levels  is  denoted  as  an  nm  factorial;  for  n  a  prime  power 
the  main  effects  and  interaction  effects  in  an  nm  factorial  are  in  a  1:1  cor¬ 
respondence  with  the  points  of  the  finite  projective  geometry  PG(m-l,n)  .  For 
2 

example,  the  n  factorial  consists  of  two  main  effects,  say  A  and  B  with 

levels  (A)  and  (B)  respectively,  i,j  =  0,1, 2, . . .  ,n-l ,  and  n-1  two  factor 
u  J  u 

interactions  AB  S,  s  =l,2,,..,n-l  with  levels  (AB  S)  +u  u  for  ui  +  ugu  » 

Ui  Us  j 

uQ,  Uj,  u2,  ... ,  u  where  the  are  elements  of  the  Galois  field  GF(n)  , 

and  the  n  +  1  effects  are  in  a  1:1  correspondence  with  the  points  of  PG(l,n)  . 

Each  of  the  n  +  1  effects  is  associated  with  a  set  of  n-1  single-degree-of-  t 

2 

freedom-contrast  parameters  making  a  total  of  (n+l)(n-l)=n  -1  parameters; 

2 

if  the  mean  is  adjoined  to  the  set  of  contrasts  then  the  n  single-degree-of- 

freedom-contrast  parameters  are  in  a  1:1  correspondence  with  the  points  of  the 

2 

finite  Euclidean  geometry  EG(2,n)  .  Therefore,  the  n  combinations  are 

2 

in  a  1:1  correspondence  with  the  n  single-degree-of-freedom-contrast  param¬ 
eters  in  EG(2,n)  . 

For  n  =  4,  the  levels  of  the  main  effects  and  interactions  are  given  by 
u  2 

(A)i ,  (B)j,  and  (AB  +y  y  ,  where  uQ  =0,^*1,  u2=x,  u3=l  +  x=x  are 

i  s  j  , 

the  marks  of  GF(4),  i,j  =  0,  1,  2,  3,  and  s  =  1,  2,  3  .  Let  (A),  be  the  rows 

and  (B)j  be  the  columns  of  a  latin  square  of  order  4  as  follows: 
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column  1  =  (B)q 

column  2  =  (B)^ 

column  3  =  (B) 

column  4  =  (B)^ 

o 

z 

I— 

II 

> 

o 

00 

01 

02 

03 

row  2  =  (A)^ 

10 

1 1 

12 

1  3 

row  3  =  (A)^ 

20 

21 

22 

23 

row  4  =  (A)3 

30 

31 

32 

33 

In  the  above  only  the  subscript  of  the  combination  u^Uj  an<^  of  the  effects  A  and 

B  is  given  for  each  row-column  intersection.  Thus,  (A)  =  (A)  consists  of 

U0 

the  n  =  4  subscripts  00,01,02,03  of  the  combinations  U0U0»  uoul»  uoU2’  U0U3  ' 

The  remaining  levels  are  similarly  defined. 

A  symbol  in  a  la  tin  square  corresponds  to  those  combinations  u.u.  for 

u 

which  u.  +  u  u  for  interaction  effect  AB  S ,  is  a  constant,  with  each  constant 
1  s  j 

corresponding  to  one  of  the  n  symbols  in  the  latin  square  of  order  n  ,  Also, 
n-1  latin  squares  of  order  n  may  be  formed  for  s  =  1,  2,  ...  ,  and  n-1;  this 
set  of  latin  squares  forms  an  0(n,n-l)  set.  For  n  =  4  the  0(4,3)  set  is 
formed  as  follows  (additional  detail  may  be  found  in  Mann  [1949] ,  chapter  VIII, 
Kempthorne  [1952],  pages  331-340,  and  Federer  [1955],  chapters  VII,  IX  and  XV): 


1 


1 

(AB  ) 


Wj 


‘2 

‘3 


00  +  1  1  +  22  +  33  -1 

01  +  10  +  23  +  32  -  II 

02  +  13  +  20  +  31  -  III 

03  +  12  +  21  +  30  “*■  IV 


00=1 

01  =  11 

02=111 

03=  IV 

10=11 

11  =  1 

12=  IV 

13=111 

20=111 

21  =  IV 

22=1 

23=11 

30=  IV 

31  =  111 

32=11 

33=1 

* 
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^u.+u 

1 


2  J 


u 


(AB 


3 


) 


U.+U,U. 

i  3  j 


ruQ  00  +  13  +  21+32 
Uj  03  +  10  +  22  +  31 
u2  01  +  12  +  20  +  33 
^  u3  02  +  1 1  +  23  +  30 


u 


1 

2 

3 


00  +  12  +  23  +  31  — 
02  +  10  +  21  +  33  - 
03  +  11  +  20  +  32  - 
01  +  13  +  22  +  30  - 


a 

a 

V 

6 

P 

P 

P 

6 

V 

a 

V 

V 

a 

P 

6 

6 

6 

P 

a 

V 

W 

w 

z 

X 

Y 

X 

X 

Y 

w 

z 

Y 

Y 

X 

z 

w 

Z 

z 

w 

Y 

X 

where  the  first  column  to  the  right  of  the  brace  represents  the  ui  obtained  from 


the  subscript. 


In  the  above  the  complete  confounding  scheme  of  sources  of  variation  in 
the  0(4,3)  set  and  the  effects  in  the  factorial  may  be  illustrated  in  the  following 


analysis  of  variance  table  wherein  the  total  sum  of  squares  has  been  orthogonally 


decomposed  into  the  sums  of  squares  related  to  the  above  confounding  scheme  as 


follows  ; 


Source  of  variation 
Correction  for  mean 


Rows  =  A  effect 


Columns  =  B  effect 


I, 


Roman  numbers  =  (AB  )  effect 

u2 

Greek  letters  =  (AB  )  effect 
u3 

Latin  letters  =  (AB  )  effect 


Degrees  of  freedom 
1 

3 

3 

3 

3 

3 


Total 


16 


Instead  of  relating  the  mutually  orthogonal  latin  squares  of  order  4  to  a 
2  4 

4  factorial  we  may  relate  them  to  a  2  factorial  in  the  following  manner,  i.  e. , 
we  consider  EG(4, 2)  and  GF(2)  with  elements  0  and  1  .  Let  the  16  row-column 

intersections  be  numbered  as  follows: 
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column 

row 

1 

2 

3 

4 

1 

0000 

0001 

0010 

0011 

2 

0100 

0101 

0110 

0111 

3 

1000 

1001 

1010 

101  1 

4 

1 100 

1101 

1 110 

11 11 

where  the  subscripts  in  the  above  table  represent  the  combination  a  b.c.d,  of 

g  h  1  j 

the  factors  a,  b,  c,  and  d  with  two  levels  (0  and  1)  each.  The  rows  correspond 
to  factorial  effects  A,  B,  and  AB  and  the  columns  correspond  to  factorial  effects 
C,  D,  and  CD  .  (This  form  of  constructing  latin  squares  has  been  used  by  Fisher 
and  Yates  [1957]  for  latin  squares  of  order  8  and  by  Federer  [1955]).  Then,  let 


the  symbols  in  the  3  latin  squares  be  represented  by  the  following  scheme: 
Factorial  generators  Combinations  latin  squares 


(AC)0,  (BD)q,  (ABCD)q 
(AC)0,  (BD)lf  (ABCDJj 
(AOj,  (BD)q,  (ABCDJj 

(AC)  j ,  (BD)1}  (ABCD)q 

(AD) q,  (ABC)0,  (BCD)q 
(AD)q,  (ABC)1}  (BCD)  j 
(AD)j ,  (ABC)0,  (BCDJj 
(AD)j ,  (ABOj,  (BCD)  Q 


0000  +  0101  +  1010  +  1  111  =  1 
0001  +  0100  +  1011  +  1110  =  II 

0010  +  0111  +  1000  +  1101  =  III 

0011  +  0110  +  1001  +  1100  =  IV 

0000  +  0110  +  101  1  +  1  101  =  w 

0010  +  0100  +  1001  +  1 1 1 1  =  X 

0001  +  0111  +  1010  +  1100  =  z 

0101  +  001  1  +  1000  +  1110  =  Y 


I 

II 

III 

IV 

II 

I 

IV 

III 

III 

IV 

I 

II 

IV 

III 

II 

I 

w 

z 

X 

Y 

X 

Y 

w 

z 

Y 

X 

z 

w 

z 

w 

Y 

X 

(ACD)0,  (BC)0,  (ABD)0 
(ACD)q,  (BOj,  (ABD)j 
(ACDJj,  (BC)0,  (ABD)l 
(ACDJj,  (BC)j  ,  (ABD)q 


0000  +  01  11  +  1110  +  1001  =  a 

1010  +  0100  +  0011  +  1101  =  (3 

1000  +  0110  +  1 1 1  1  +  0001  =  y 

0010  +  0101  +  1011  +  1100  =  6 


a 

V 

6 

P 

P 

5 

V 

a 

V 

a 

P 

6 

6 

a 

V 

*  Note:  Some  authors  use  lower  case  letters  to  denote  the  factors  and  capital 
letters  to  denote  effects  or  levels  of  effects;  we  follow  that  usage  here. 
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The  correspondence  of  the  latin  squares  obtained  from  complete  confound- 

2  4 

ing  considering  a  4  factorial  and  considering  a  2 

factorial  is  demonstrated 

in  the  following  analysis  of  variance  table: 

Source  of  variation 

dearees  of  freedom 

Correction  for  mean 

1 

Rows  =  A  effect  in  r  factorial 

3 

A  effect  in  2  4  factorial 

1 

B  "  "  24  " 

1 

AB  "  "  24 

1 

2 

Columns  =  B  effect  in  4  factorial 

3 

C  effect  in  24  factorial 

1 

D  "  "  24  " 

1 

CD  "  "  2 4  " 

1 

v-i  2 

Roman  numbers  r-  effect  in  4  factorial 

3 

AC  effect  in  2 4  factorial 

1 

BD  "  "  24  " 

1 

ABCD  "  "  24  " 

1 

u2  ? 

Greek  letters  =  AB  effect  in  4C  factorial 

3 

t 

ACD  effect  in  24  factorial 

1 

BC  "  "  24  " 

1 

ABD  "  « 24 

1 

u3  2 

Latin  letters  =  AB  effect  in  4  factorial 

3 

AD  effect  in  24  factorial 

1 

ABC  "  " 24 

1 

BCD  "  "  24  " 

1 

Total  16 


4 

It  should  be  noted  here  that  the  effects  in  the  2  map  directly  into  the 
2  2 

4  projective  geometry  or  PG(1,2  ).  Likewise,  even  though  one  more  set  of 

* 

generators  is  available,  viz. 

t 
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Generators 

interaction 

Roman  numbers 

=  AD,  BC 

ABCD 

Greek  letters 

=  AC,  ABD 

BCD 

Latin  letters 

=  BD,  ABC 

ACD 

the  three  orthogonal  latin  squares  produced  are  the  same  ones.  Since  the  third 

effect  above  is  obtained  as  the  product  of  two  generators  (exponents  mod  2)  we 

need  consider  only  two  generators.  Multiplying  these  by  CD  (exponents  mod  2) 

we  obtain  the  generators  of  the  preceding  scheme.  Hence,  even  though  two  dif  - 

ferent  complete  confounding  schemes  are  available  there  is  a  simple  one-to-one 

mapping  of  one  set  into  the  other  set.  Although  nothing  interesting  turns  up  here, 

it  would  be  interesting  to  study  the  various  complete  confounding  schemes  in  the 

4 

latin  square  of  order  9  as  related  to  the  3  factorial. 

As  a  second  illustration  of  the  use  of  complete  confounding  to  construct 

latin  squares,  let  us  consider  a  latin  square  of  order  6  .  Using  the  notation 
and  concepts  of  Raktoe  [1969]  on  mixed  prime  factorials  as  related  to  rings  and 

2  2  2 

elements  of  ideals  in  the  rings  we  designate  the  6  as  a  2  (3)  factorial  and 

represent  a  combination  by  ghij  where  g,  h  are  members  of  the  ideal  I  (3)  and 

2  2 

i,  j  are  members  of  the  ideal  I  (4)  .  The  effects  in  the  2  and  in  the  3  factorials 
are  denoted  respectively  by: 


4 


D' 


A3  3 
A  B 


4  4 
C  D 


4  2 
C  D 


The  remaining  interactions  are  given  below  in  the  analysis  of  variance  table: 
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Source  of  variation 


Degrees  of  freedom 


Correction  for  mean 

1 

3  4 

Rows  =  A  C 

5 

A3 

c4 

A3  X  C4 

„  ,  3  4 

Columns  =  B  D 

5 

B3 

D4 

B3  X  D4 

3  3*4 

Treatments  or  symbols  =  A  B  C  D 

5 

3  3 

A  B 

4  4 

C  D 

.  3_3  4  4 

A  B  X  C  D 

Remainder 

20 

c4d2 

2 

A3  X  D4 

2 

a3x  c4d4 

2 

A3  X  C4D2 

2 

B3  X  C4 

2 

B3  X  C4D4 

2 

B3  X  C4D2 

2 

3  3  4 

A  B  X  C 

2 

.  3_  3  4 

A  B  X  D 

2 

.  3_3  _4  2 

A  B  X  C  D 

2 

Total 


36 


Let  us  now  set  up  the  6  rows  and  the  6  columns  of  a  latin  square  of  order 
6  with  the  corresponding  designation  of  the  36  combinations  as  follows: 

Columns 


Rows 

(b3d4)0 

(b3d4)1 

<b3d4)2 

(B3D4)3 

(bV)4 

(bV,5 

<a3c\ 

0000 

0304 

0002 

0300 

0004 

0302 

(A3c4)1 

3040 

3344 

3042 

3340 

3044 

3342 

(a3c4)2 

0020 

0324 

0022 

0320 

0024 

0322 

<a3c4)3 

3000 

3304 

3002 

3300 

3004 

3302 

(A3c4)4 

0040 

0344 

0042 

0340 

0044 

0342 

(a3c4)5 

3020 

3324 

3022 

3320 

3024 

3322 

Now  let  the  levels  of  A  B  C  D  correspond  to  the  symbols  in  a  latin 
square  of  order  6  as  follows: 


Levels  Combination  for  which  3q+3h+4i-t4j .  mod  6.  is  constant  Symbol 


(A3B3C4D4)q  0000  +  3342  +  0024  +  3300  +  0042  +  3324  -  0 

3  3  4  4 

(A  B  C  D  )  0304  +  3040  +  0322  +  3004  +  0340  +  3022  -*  1 

3  3  4  4 

(A  B  C  D  )  0002  •!•  3344  +  0020  +  3302  +  0044  +  3320  -  2 

3  3  4  4 

(A  B  C  D  )  0300  +  3042  +  0324  +  3000  +  0342  +  3024  -  3 

(A3B3C4D4)4  0004  +  3340  +  0022  +  3304  +  0040  +  3322  -  4 

3  3  4  4 

(A  B  C  D  L  0302  +  3044  +  0320  +  3002  +  0344  +  3020  -  5 


#1030 


-11- 


This  produces  the  following  latin  square  of  order  6 


0 

1 

2 

3 

D 

5 

1 

2 

3 

D 

5 

0 

2 

3 

D 

5 

0 

1 

3 

D 

5 

0 

1 

2 

D 

5 

0 

1 

2 

3 

5 

0 

1 

■  1  • 

2 

3 

m 

Alternatively  we  could  have  used  levels  of 


.3  3  4  2 
A  B  C  D 


to  construct  the 


following  latin  square  of  order  6: 


Levels 


(a3b3c4d2) 

(aVcV) 

..3  3  4  2, 
(A  B  C  D  ). 

..3  3  4  2/ 
(A  B  C  D  ). 

/R3  3  4  2  ' 
(A  B  C  D  )i 

..3-3  _^4  2 . 
(A  B  C  D  ). 


1 


Combinations  for  which  3g+3h+41+2j ,  mod  fe.is  constant  Symbol 


0000  +  3344  +  0022  +  3300  +  0044  +  3322  0 

0302  +  3040  +  0324  +  3002  +  0340  +  3024  1 

0004  +  3342  +  0020  +  3304  +  0042  +  3320  2 

0300  +  3044  +  0322  +  3000  +  0344  +  3022  3 

0002  +  3340  +  0024  +  3302  +  0040  +  3324  4 

0304  +  3042  +  0320  +  3004  +  0342  +  3020  5 


latin  square  of  order  6 


0 

5 

4 

3 

2 

1 

1 

0 

5 

4 

3 

2 

2 

1 

0 

5 

4 

3 

3 

2 

1 

0 

5 

4 

4 

3 

2 

1 

0 

5 

5 

4 

3 

2 

1 

0 
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Thus,  the  above  square  is  simply  a  column  permutation  of  the  previous  one.  As 
there  are  no  other  sets  of  5  degrees  of  freedom  leading  to  a  latin  square  of  order 

6  (i.e.  A3,  B3,  and  A3B3  exhaust  the  three  single  degrees  of  freedom  from  the 

2  4  4  4  4  4  2 

2  factorial  and  C  ,  D  ,  C  D  ,  and  C  D  exhaust  all  sets  of  2  degrees  of 

2 

freedom  from  the  3  factorial,),  it  is  not  possible  to  obtain  a  latin  square  of 
order  6  orthogonal  to  either  of  the  preceding  ones  using  complete  confounding 
schemes. 

For  a  latin  square  of  order  10  we  may  use  levels  of  ABCD,ABCD, 
A5B5C6D8,  or  A5B5C6D4  to  form  four  different  latin  squares  of  order  10. 

II.  2.  Partial  Confounding 

In  the  last  section  use  was  made  of  complete  confounding  of  effects  in  a 
factorial  with  the  rows,  columns,  and  symbols  in  a  latin  square.  In  this  section 
some  of  the  factorial  effects  will  be  partially  confounded  with  row  (column  or 
symbol)  contrasts,  i.e.  contrasts  among  levels  of  an  effect  will  be  completely 
confounded  with  a  subset  of  the  row  (column  or  symbol)  contrasts  and  will  be 
unconfounded  with  the  remaining  contrasts,  and  vice  versa.  In  complete  con¬ 
founding  no  subset  of  contrasts  among  the  levels  of  a  factorial  effect  can  be 
separated  from  contrasts  among  the  rows  (columns  or  symbols).  (See,  e.g., 

Yates  [1937]  and  Federer  [1955]  ).  For  example,  the  latin  square  of  order  4 

4 

could  be  considered  as  a  2  factorial  as  in  the  preceding  section,  with  the 
following  scheme  of  confounding: 

Columns 


Rows 

ii 

Q 

o 

2=  (C)j 

u> 

II 

8 

o 

II 

§ 

1 

(A)0,(B)0 

0000 

0011 

0010 

0001 

2 

0101 

0110 

0100 

0111 

3 

(A)i,(B)0 

1000 

1011 

1010 

1001 

4 

(Alj^Blj 

1101 

1 110 

1100 

1 1 1 1 
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If  we  set  up  the  latin  square  symbols  for  the  above  as 


then 


a 

P 

V 

6 

P 

a 

6 

V 

\ 

6 

P 

a 

6 

a 

P 

the  symbols  correspond  to  the  following  combinations: 

at  0000  +  01  10  +  1001  +  1100  =  (ABCD)Q  +  other  effects 

(3:  0011  +  0101  +  1010+1111  =  (ABCD)Q  + 

y:  1000  +  1  1  10  +  0010  +  011  1  =  (ABCD)  + 

6:  0001  +  0100  +  101  1  +  1  101  =  (ABCD)  + 


It  is  known  that  this  latin  square  has  no  ^r*hogonal  mate  (Hedayat  [1969]  ). 
This  means  that  no  orthogonal  partition  of  the  remaining  sum  of  squares  can  be 
made  which  forms  a  latin  square. 


If  on  the  other  hand,  the  latin  square  used  is 


a 

P 

V 

6 

P 

a 

5 

V 

V 

6 

a 

P 

6 

P 

a 

the  combinations 


corresponding  to  the  Greek  letters  are: 


a:  0000  +  0110  +  1010+  1111  =  (ABCD)Q  +  other  effects 

(3:  001  1  +  0101  +  1001  +  1  100  =  (ABCD)Q  + 

y:  0010  +  01  11  +  1000  +  1  1  10  =  (ABCD) {  +  (AC^  +  other  effects 

6:  0001  +  0100  +  101  1  +  1  1 0 1  =  (ABCD)  {  +  other  effects 
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This  square  has  two  mutually  orthogonal  mates  and  hence  there  must  be  partitions 
of  the  sums  of  squares  into  orthogonal  components  which  correspond  to  the  symbols 
in  a  latin  square. 

Instead  of  inserting  symbols  in  the  latin  square  of  order  4,  denote  the 

symbols  in  the  latin  square  by  the  following  partial  confounding  scheme,  where 

a  fractional  replicate  is  a  subset  of  a  complete  factorial: 

i)  add  the  two  1/8  replicates  generated  by  ((A)Q,  (D)Q,  (BC)Q)  and 

((A)  ,  (C)  ,  (ABD)j)  to  obtain  the  4  combinations  (0000  +  01  10)  + 

(1010  +  1111)  and  denote  these  4  combinations  as  symbol  a  , 

ii)  add  the  two  1/8  replicates  generated  by  ((D)^,  (AB)^,  (AC)Q)  and 

((AB)q,  (C)q,  (AD)j)  to  obtain  combinations  (0101  +  1011)  +  (1100  + 

0001)  and  denote  these  4  combinations  as  symbol  (3  , 

iii)  add  the  two  1/8  replicates  generated  by  ((A)p  (D)Q,  (ABC)^)  and 
((A)q,  (C)j,  (BD)0)  to  obtain  combinations  (1000  +  1110)  +  (0010  + 

0111)  and  denote  these  4  as  symbol  v  > 

iv)  add  the  two  1/8  replicates  generated  by  ((AB)Q,  (AC)^,  (D)^)  and 
((AB)p  (C)Q ,  (BD)j)  to  obtain  the  combinations  (1101  +  0011)  + 

(0100  +  1001)  and  denote  these  4  as  symbol  6  . 

This  procedure  results  in  the  following  latin  square  of  order  4: 


a 

6 

Y 

|3 

P 

a 

6 

Y 

V 

P 

a 

6 

6 

Y 

(3 

a 
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Obviously,  one  could  take  any  pair  of  1/8  replicates  such  that  the  4  combina¬ 
tions  are  in  different  rows  and  in  different  columns  to  form  the  combinations  for 
a  given  symbol. 

The  above  type  of  partial  confounding  results  in  the  class  of  latin  squares 
denoted  as  half-plaid  latin  squares  (See  Federer  [1955]  chapters  IX  and  XV  and 
Yates  [1937] ).  If  partial  confounding  were  utilized  in  rows  as  well  as  in  columns 
the  resulting  square  would  be  denoted  as  a  plaid  latin  square  (so-called  because 
of  its  resemblence  to  plaid  cloth  if  the  effects  confounded  were  of  different  colors). 
The  three  types  of  squares  are  illustrated  below  for  a  latin  square  of  order  6  where 
the  factorial  effects  are  as  described  in  statistics  books  (e.  g.  ,  Federer  [1955]): 


Complete  confounding  of  effects 
Columns 


1  = 

2  = 

3  = 

4  = 

5  = 

6  = 

Rows 

(A»0’,C)0 

(A)0,(O2 

WpICto 

(A)p  (C)x 

(A)2,(o2 

1  =  (B)0’(D)0 

0000 

0010 

0020 

1000 

1010 

1020 

2  =  (B)0,(D)1 

0001 

0011 

0021 

1001 

1011 

1021 

3=  (B)0,(D)2 

0002 

0012 

0022 

1002 

1012 

1022 

4  =  (B)  j ,  (D)q 

0100 

0110 

0120 

1  100 

1 110 

1120 

5  =  (B)1,(D)1 

0101 

0111 

0121 

1  101 

1  111 

1121 

6  =  (B)t ,  (D)2 

0102 

0112 

0122 

1102 

1112 

1122 

Partial  confounding  of  effects  with  columns 
Columns 


Rows 

1  *  <c>0 

2  =  (C)j 

3  =  (C)2 

4  =  (CD)q 

5  =  (CD)  j 

6  =  (CD)2 

1  -  (B)o,(D)0 

0000 

0010 

0020 

1000 

1010 

1020 

2  = 

0001 

0011 

0021 

1021 

106l 

mmm 

3  =(B)0,(D)2 

0002 

0012 

0022 

1012 

■ 

4  =  (B^,^ 

1  100 

1 1  10 

1120 

0100 

0110 

0120 

5  =  (B)1,(D)1 

1 101 

1 1  1 1 

1121 

0121 

0101 

0111 

6  = 

1 102 

1  1  12 

1122 

0112 

0122 

0102 
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Partial  confounding  in  both  rows  and  columns 


Columns 


Rows 

M 

o 

o 

2  =  (C)1 

3  =  (C)2 

4  =  (CD)0 

5  =  (CD)J 

6  =  (CD)2 

ii 

o 

o 

00 

10 

20 

00 

10 

20 

2  =  (D)] 

01 

11 

21 

21 

01 

1  1 

3  =  (D)2 

02 

12 

22 

12 

22 

02 

4  =  (CD2)q 

00 

11 

22 

00 

22 

1  1 

5  =  (CD*)l 

02 

10 

21 

- 

21 

10 

02 

MESiM 

01 

12 

20 

12 

01 

20 

In  the  last  table  above  only  the  subscripts  for  combinat'ans  of  factors  c 
and  d  have  been  inserted.  There  is  some  difficulty  in  inserting  subscripts  for 
factors  a  and  b  such  that  these  effects  are  orthogonal  to  both  rows  and  columns. 
In  any  event,  this  problem  requires  further  study  to  determine  if  half-plaid  latin 
squares  and  plaid  latin  squares  lead  to  latin  squares  not  of  the  same  type  as  given 
by  complete  confounding.  If  the  three  types  of  latin  squares  of  order  6  can  be 
produced  by  partial  and  complete  confounding,  this  would  be  an  interesting  result. 
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III.  Fractional  Replication  Construction  of  Q(n,  t)  Sets 


Any  latin  square  may  be  considered  as  an  n  fraction  of  an  n  factorial 
where  the  rows  represent  levels  of  one  factor,  the  columns  represent  the  levels 
of  the  second  factor,  and  the  symbols  in  the  latin  square  represent  the  levels  of 
the  third  factor.  As  an  Illustration,  consider  the  latin  square  of  order  3  where 
the  9  combinations  represent  the  1/3  fraction  of  a  3  factorial  as  follows: 


Columns 


Rows 

0 

1 

2 

0 

000 

012 

021 

1 

102 

111 

120 

2 

201 

210 

222 

3 

The  above  is  the  1/  3  fraction  of  a  3  corresponding  to  (ABC1  .  .  „  ,  , 

h+i+j  =  0,mod  3  . 

Since  this  is  a  regular  fraction  we  may  write  out  the  aliasing  structure  in  this 

i 

fraction  as  follows; 


M  +  ABC 

A  +  AB2C2  +  BC 

B  +  AB2C  +  AC 
2 

C  +  ABC  +•  AB 

2  2  2 
AB  +  AC  +  SC 


where  the  effects  connected  with  a  plus  sign  are  completely  confounded  with  each 

other.  In  the  above  latin  square  the  symbols  0,1,2  correspond  to  the  levels  of 

the  third  factor,  c.  Now  if  we  set  up  a  second  latin  square  in  which  the  symbols, 

2 

say  a? ,  (3 ,  \  ,  correspond  to  the  levels  of  AB  ,  the  resulting  square  will  be 

2 

orthogonal  to  the  first  one.  The  square  corresponding  to  levels  of  (AB  ),  , .  ,  , 

i+2j ,  mod  3 


a 

V 

(3 

(3 

a 

V 

\ 

P 

a 

) 
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000  +  1  1  1  +  222  =  a 
021  +  210  +  102  =  p 
201  +  012  +  120  =  y 
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The  class  of  fractional  replicates  constituted  as  an  n  fraction  of  an  n 

factorial  becomes  an  important  one  to  study  as  it  relates  to  construction  of  mutu- 

-3  9 

ally  orthogonal  latin  squares.  In  particular,  all  2  fractions  of  a  2  and  all 
-2  6 

3  fractions  of  a  3  with  all  possible  aliasing  structures  could  produce  several 
sets  of  mutually  orthogonal  latin  squares.  This  could  have  interesting  conse¬ 
quences  in  finite  geometry. 

The  structure  of  the  left-hand  set  of  parameters  in  an  aliasing  structure 
will  have  a  pattern;  for  example,  for  n  =  4,  5,  and  7,  the  patterns  are: 

n  =  4  n  =  5  n  =  7 

M  +  ABC  M  +  ABC  M  +  ABC 


A 

A 

A 

B 

B 

B 

C 

C 

C 

AB2 

AB2 

AB2 

AB3 

AB3 

AB3 

AB4 

AB4 

AB5 

AB6 

Note  that  although  ABC  was  completely  confounded  with  the  mean,  any  one  of 

U  V 

the  other  three-factor  interaction  components  AB  C  ,  u,v  =  l,2,...,n-l  could 

have  been  utilized  equally  well.  Also,  note  that  the  levels  of  C  corresponding 

to  symbols  produce  a  latin  square  ,  and  that  the  levels  of  effects  below  the 

factor  B  produce  a  set  of  n-1  mutually  orthogonal  latin  squares. 

-1  3 

In  general  we  want  to  look  at  all  possible  n  fractions  of  an  n  factorial, 
i3> 


i.  e. ,  the  subset  of 


0 


combinations  for  which  the  levels  of  C  are  the  symbols 


in  a  latin  square  and  to  study  their  patterns  especially  for  n  =  7,8,  and  9.  All 


possible  fractions,  or  rather  all  forms  of  the  aliasing  structure,  could  be  classi' 
fied  into  all  types  of  t  mutually  orthogonal  latin  squares,  0(n,t)  for  t  =  1,2, 


. . .  ,n-l  .  Perhaps  this  is  the  manner  in  which  the  geometries  of  various  values 
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of  n  can  be  exhaustively  studied.  In  fractional  factorial  notation  we  want  to 
study  all  possible  aliasing  patterns  for  one  latin  square,  for  two  latin  squares, 
etc.  as  given  by: 


etc. 


2 

where  J3Q  is  the  n  +  n-2  vector  containing  the  interaction  effect  parameters, 

W  is  the  4  X  (n2  +  n  -  2)  matrix  of  aliasing  coefficients,  X  is  one  of  the 

two  factor  interaction  effects  in  corresponding  to  a  column  of  zero  co- 

*  * 

efficients  in  W,  and  and  W  correspond  to  g,Q  and  W  with  the  param¬ 
eter  X  deleted.  For  n  =  3,  ^  =  (AB,  AB2,  AC,  AC2,  BC,  BC2,  ABC,  ABC2  , 
AB2C,  AB2C2)  and  W  is  equal  to 

/0  0  0  0  0  0  1  0  0  o\ 
0000100001 
0010000010 
\l  000000  1  00/ 

2 

Since  there  are  three  columns  containing  all  zeros  X  could  be  either  AB  , 

2  2  2  * 

AC  ,  or  BC  .  Selecting  X  as  AB  ,  say,  there  would  be  no  columns  in  W 

which  contain  all  zeros.  Thus,  to  obtain  an  0(n,  2)  set  from  a  given  0(n,  1) 

set,  at  least  one  column  in  W  should  be  all  zeros .  Likewise,  to  obtain  an 

$ 

0(n,3)  set  from  a  given  0(n,  2)  set,  at  least  one  column  in  W  should  contain 
all  zeros. 
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We  now  wish  to  illustrate  the  use  of  fractional  replication  procedures  to 
construct  latin  squares  which  are  mateless  and  which  have  orthogonal  mates. 


To  illustrate  let  us  consider  the  four  standard  latin  squares  of  order  4  which  are 
(Fisher  and  Yates  [1957]  ): 


0 

1 

2 

3 

1 

2 

3 

0 

2 

3 

0 

1 

3 

0 

1 

2 

Square  I 


0 

1 

2 

3 

1 

3 

0 

2 

2 

0 

3 

1 

3 

2 

1 

0 

Square  II 


0 

1 

2 

3 

1 

0 

3 

2 

2 

3 

1 

0 

3 

2 

0 

1 

Square  III 


0 

1 

2 

3 

1 

0 

3 

2 

2 

3 

0 

1 

3 

2 

1 

0 

Square  IV 


It  is  known  (Hedayat  [1969]  )  that  the  first  three  squares  are  mateless  (There  is 
no  transversal  through  000. )  and  that  the  last  square  belongs  to  an  0(4,  3)  set. 

Now  number  the  rows  as  0, 1,  2,  3  and  denote  these  as  levels  of  the  factor 
a;  number  the  columns  as  0,1,  2,3,  and  denote  these  as  levels  of  factor  b, 
and  the  symbols  in  the  latin  squares  by  0,1,  2,  3,  the  levels  of  the  factor  c  . 
Then,  in  factorial  notation  the  above  16  combinations  form  a  one-fourth  fraction 

3 

of  a  4  factorial  treatment  design.  The  aliasing  scheme  for  the  fractional 

replicate  given  as  square  IV  is 

»  M  +  ABC 

A  +  BC  +  AB2C2  +  AB3C3 
B  +  AC  +  AB2C  +  AB3C 
C  +AB  +  ABC2  +  ABC3 


where  u^  =  1,  2  means  u2,  and  3  means  u3  from  GF(4)  and  where  the 

effects  connected  with  a  plus  sign  are  completely  confounded  with  each  other. 
The  completion  of  the  remaining  two  aliasing  structures  results  in  the  complete 

-1  3 

aliasing  structures  for  this  4  fraction  of  the  4  factorial;  these  two  are; 
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2  3  2  3  2 

AB  +  AC  +  BC  +  AB  C 

3  2  3  2  3 

AB  +  AC  +  BC  +  AB  C 

2  3 

If  we  use  the  levels  of  AB  and  of  AB  to  form  two  latin  squares,  these  two 
with  square  IV  form  an  0(4,3)  set  of  mutually  orthogonal  latin  squares. 

Now,  let  us  return  to  the  set  of  four  standard  squares  given  above  and 
we  note  that  only  four  combinations  in  square  IV  are  replaced  to  obtain  squares 
I,  II,  and  HI.  These  are: 


additional  combinations 

combinations  replaced  in  IV 

Square  I 

112,  130,  310,  332 

110,  132,  312,  330 

"  II 

113,  120,  210,  223 

110,  123,  213,  220 

"  III 

213,  230,  320,  331 

220,  231,  321,  330 

The  aliasing  structure  (without  the  coefficients)  is  given  on  the  following  page  for 
all  four  standard  latin  squares  of  order  4.  The  l/4  replicate  given  by  square  IV 
forms  a  regular  fraction.  The  remaining  three  fractional  replicates  are  such  that 
none  of  the  additional  effects  are  unconfounded  with  the  effects  M,  A,  B,  or  C  of 
the  original  latin  squares  of  order  4.  Since  this  is  true  no  linear  combination 
of  these  effects  will  be  unconfounded.  In  order  to  form  a  latin  square  which  is 
orthogonal  to  the  given  one  it  is  necessary  that  there  be  a  set  of  effects  which 
is  unconfounded  with  the  effects  in  the  given  square.  This  is  impossible  for  the 
three  squares  I,  II,  and  III  and  hence  the  squares  are  mateless,  as  is  well-known. 

It  would  be  interesting  to  ascertain  the  aliasing  structures  for  the  six 
standard  latin  squares  of  order  5  belonging  to  the  0(5,4)  set  and  for  the  fifty 
standard  latin  squares  of  order  5  for  which  are  known  to  be  mateless  (Hedayat 
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AB2  +  AC3  +  BC2  +  AB3C2 

3  2  3  2  3 

AB  +  AC  +  BC  +  AB  C 

2  3 

If  wo  use  the  levels  of  AB  and  of  AB  to  form  two  latin  squares,  these  two 
with  square  IV  form  an  0(4,  3)  set  of  mutually  orthogonal  latin  squares. 

Now,  let  us  return  to  the  set  of  four  standard  squares  given  above  and 
we  note  that  only  four  combinations  in  square  IV  are  replaced  to  obtain  squares 
I,  II,  and  III.  These  are: 


additional  combinations 

combinations  replaced  in  IV 

Square  I 

112,  130,  310,  332 

110,  132,  312,  330 

"  II 

113,  120,  210,  223 

110,  123,  213,  220 

"  III 

213,  230,  320,  331 

220,  231,  321,  330 

The  aliasing  structure  (without  the  coefficients)  is  given  on  the  following  page  for 
all  four  standard  latin  squares  of  order  4.  The  1/ 4  replicate  given  by  square  IV 
forms  a  regular  fraction.  The  remaining  three  fractional  replicates  are  such  that 
none  of  the  additional  effects  are  unconfounded  with  the  effects  M,  A,  B,  or  C  of 
the  original  latin  squares  of  order  4.  Since  this  is  true  no  linear  combination 
of  these  effects  will  be  unconfounded.  In  order  to  form  a  latin  square  which  is 
orthogonal  to  the  given  one  it  is  necessary  that  there  be  a  set  of  effects  which 
is  unconfounded  with  the  effects  in  the  given  square.  This  is  impossible  for  the 
three  squares  I,  II,  and  III  and  hence  the  squares  are  mateless,  as  is  well-known. 

It  would  be  interesting  to  ascertain  the  aliasing  structures  for  the  six 
standard  latin  squares  of  order  5  belonging  to  the  0(5,4)  set  and  for  the  fifty 
standard  latin  squares  of  order  5  for  which  are  known  to  be  mateless  (Hedayat 
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[  19 69 J  ).  After  a  study  of  these  fractions,  one  should  continue  such  a  study  for 

6*2 

n  =  7,  8,  and  9.  It  is  suggested  that  one  consider  a  2  fraction  instead  of  a 

3-1  9—3  3—1 

4  fraction  for  n  =  4  and  a  2  fraction  instead  of  an  8  fraction  for 

n  =  8  .  The  reason  for  this  is  that  there  is  much  more  theory  available  for 

s  =  2  in  the  sm  series  than  for  any  other  value  of  s  .  Also,  one  may  use  the 

generalized  defining  contrast  which  has  been  developed  by  Raktoe  and  Federer 

[1969]  to  a  considerable  advantage  in  writing  out  aliasing  structures  in  these 

cases.  Investigation  of  the  regular  and  irregular  fractional  replicates  obtainable 

for  various  values  of  n  could  lead  to  considerable  advances  in  the  theory  of 

mutually  orthogonal  latin  squares. 
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IV.  ANOVA  Construction  of  Q(n,t)  Sets 


There  should  be  some  procedure  which  would  utilize  the  orthogonality  of 
single  degree  of  freedom  contrasts  in  the  analysis  of  variance  (ANOVA)  and 
which  could  be  utilized  to  construct  orthogonal  latin  squares.  For  example,  one 
could  make  use  of  orthogonal  polynomial  coefficients  for  row  and  column  con¬ 
trasts  and  then  construct  mutually  orthogonal  latin  squares  from  these.  To  illus¬ 
trate,  consider  the  latin  square  of  order  4  used  previously  wherein  the  row-column 

4 

intersections  are  numbered  as  a  2  factorial,  i.  e. : 


Column 


Row 

1 

2 

3 

4 

1 

0000 

0001 

0010 

001  1 

2 

0100 

0101 

0110 

0111 

3 

1000 

1001 

1010 

1011 

4 

1100 

1  101 

1  110 

mi 

The  relation  between  the  16  contrasts  using  orthogonal  polynomial  coefficients 

4 

and  the  2  factorial  is  given  below,  where  R^,  R^,  and  R^,  are  linear, 
quadratic,  and  cubic  polynomial  contrasts  among  rows  and  CL>  Cq,  and 
are  linear,  quadratic,  and  cubic  polynomial  contrasts  among  the  columns: 
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Source  of  variation 
C.F.  M. 

Row  contrasts 


df 

1 

3 


u. 


Latin  letters  =  (AB  ) 


AD  -  2R  C  -2R  C  -R  C_+ 
L  L  O  O  L  U 


4rccl 


ABC  =  Rq(-Cl-2Cc) 
BCD  =  (-2Rl  +  rc)Cq 


A  =  -R  -  2R 

1 

L  C 

B  =  -2R  +  ft 

,i 

L  C 

AB  =  R_ 

1  - 

Q 

Column  contrasts 

C  =  -C-  2C 

"N 

1 

L  C 

D  =  -2Ct  +  C^ 

1 

L  C 

CD=  CQ 

U1 

Roman  numbers  =  (AB  ) 

AC  =  RTCT  +  4RC. 

■N 

1 

L  L  C  C 

BD  =  4R  CT  +  R^,C_ 

1  > 

L  L  C  C 

ABCD  =  R_C„ 

1 

Q  Q 

u2 

Greek  letters  =  (AB  ) 

ABD  =  -2RC  +  R„  C  _ 

r 

Q  L  Q  C 

BC  =  2RlCl  -  2RCCC+ 

4R  C  -R_C_ 

i 

L  C  CL 
ACD=  (-Rl-2Rc)Cq 

i  ^ 

1  Rows  linear  =  R  =  A  +  2B 

Li 

"  quadratic  =  Rp  =  AB 
"  cubic  =  R  =  2A  -  B 


r 


1  Columns  linear  =  CT  =  C  +  2D 

Li 


r. 


y 


>  \ 


quadratic  =  Cq  =  CD 

cubic  =  C  =  2C-D 
o 


rlcl 

RCCC 


rqcq 


rlcq 


rqcc 

rccl 


rlcc 

rqcl 

RC  CQ 


Total 


16 
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Combination 


The  individual  degree  of  freedom  contrast  matrix  for  the  above  16 


combination  is: 


p— t 

+ 

CO 

+ 

1 

CO 

+ 

1 

O' 

CO 

CO 

1 

ro 

+ 

i 

CO 

1 

i 

+ 

o 

zz 

+ 

fO 

+ 

1 

+ 

i 

ro 

ro 

ro 

1 

O' 

+ 

1 

CO 

i 

+ 

CO 

1 

-H 

o 

+ 

co 

1 

i 

i 

CO 

1 

co 

1 

co 

1 

O' 

1 

t 

1 

co 

1 

+ 

+ 

ro 

o 

O 

«— 4 

+ 

ro 

+ 

1 

CO 

1 

+ 

+ 

O' 

1 

co 

CO 

co 

1 

+ 

+ 

ro 

i 

1 

1011 

+ 

+ 

i 

CO 

CO 

+ 

i 

CO 

+ 

i 

co 

1 

1 

+ 

O' 

ro 

ro 

1 

1010 

+ 

+ 

i 

CO 

+ 

i 

ro 

+ 

i 

CO 

i 

+ 

CO 

1 

CO 

ro 

1 

O' 

f-H 

o 

o 

p*H 

+ 

+ 

i 

CO 

i 

i 

ro 

1 

i 

i 

co 

1 

+ 

+ 

co 

co 

ro 

1 

O' 

1 

O 

o 

o 

+ 

+ 

i 

ro 

CO 

1 

+ 

+ 

ro 

1 

+ 

+ 

ro 

1 

i 

O' 

i 

co 

ro 

0111 

4- 

i 

i 

CO 

ro 

4- 

i 

ro 

1 

■ 

4- 

m 

i 

' 

+ 

O' 

1 

ro 

1 

ro 

0110 

+ 

• 

■ 

ro 

+ 

i 

CO 

i 

+ 

co 

1 

i 

+ 

ro 

1 

ro 

1 

co 

O' 

1 

o 

O 

+ 

i 

' 

ro 

1 

1 

i 

co 

1 

+ 

+ 

co 

+ 

+ 

ro 

co 

ro 

O' 

o 

o 

o 

+ 

i 

1 

ro 

ro 

+ 

+ 

fO 

i 

1 

CO 

’ 

1 

O' 

co 

1 

CO 

0011 

+ 

CO 

+ 

+ 

CO 

+ 

i 

O' 

1 

co 

1 

co 

ro 

+ 

■ 

CO 

+ 

• 

0010 

+ 

CO 

1 

+ 

+ 

+ 

i 

co 

co 

fO 

T 

+ 

' 

ro 

+ 

i 

ro 

o 

o 

o 

+ 

CO 

+ 

+ 

' 

i 

co 

co 

co 

O' 

1 

' 

co 

1 

1 

ro 

o 

o 

o 

o 

+ 

CO 

1 

+ 

+ 

ro 

+ 

+ 

O' 

co 

1 

CO 

CO 

+ 

+ 

co 

+ 

+ 

Contrast 

Mean 

rl 

oP 

o 

a ; 

i-j 

o 

o 

o 

R  C 

L  L 

0£ 

o 

o 

l-J 

a; 

o 

a? 

a 

o 

o 

a 

o 

o 

a 

a; 

o 

o 

a; 

a 

o 

o 

a; 

o 

o 

o 

a: 
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The  corresponding  single  degree  of  freedom  contrast  matrix  for  the  2  factorial  is: 


Source  of  variation 

df 

C.F.  M. 

1 

Row  contrasts 

3 

r 

=  A  +  2B 

A  =  -R  -  2R 

L  O 

1 

1  Rows  linear  =  R 

L 

B  =  -2R  +  R 

L  O 

1 

>  < 

1  "  quadratic  = 

O 

ii 

£ 

> 

00 

II 

O 

1  "  cubic  =  R_, 

k.  C 

=  2A  -  B 

Column  contrasts 


u. 


Latin  letters  =  (AB  ) 


AD  =  2RlCl-2RcCc-RlCc+ 

4rccl 

ABC  =  Rq(-Cl-  2Cc) 

BCD  =  (-2Rl  +  RC)CQ 


C  =  -CT-  2C_ 

> 

1 

L  C 

D  =  -2Ct  +  C^ 

1  > 

L  C 

CD  -  CQ 

U1 

Roman  numbers  =  (AB  ) 

AC  =  RT  CT  +  4RC 

*\ 

1 

L  L  C  C 

BD  =  4R  CT  +  R_C_ 

1  > 

L  L  C  C 

ABCD  =  R„C_ 

1 

Q  Q 

U2 

Greek  letters  =  (AB  ) 

ABD  =  -2RC.  +  R~C 

r 

Q  L  Q  C 

BC  =  2RlCl  -  2RCCC+ 

4R  C  -R^Ct 

} 

1 

L  C  CL 
ACD=  (-Rl-2Rc)Cq 

1  > 

1  Columns  linear  =  C^  =  C  +  2D 

"  quadratic  =  Cq  =  CD 

"  cubic  =  C_  =  2C-D 

o 


rlcl 

Rccc 

Roco 


1 


■>  N- 


1 


rlcq 


rqcc 

rccl 


rlcc 

rqcl 

RC  CQ 


Total 


16 
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The  corresponding  single  degree  of  freedom  contrast  matrix  for  the  2  factorial  is: 


The  corresponding  single  degree  of  freedom  contrast  matrix  for  the  2^ 
factorial  is: 


#1030 


The  particular  contrast  matrix  utilized  is  not  unique,  as  has  been 
demonstrated  above.  All  orthogonal  contrast  matrices  resulting  in  latin  squares 
could  be  considered.  For  example,  other  sets  of  contrasts  among  rows  (or  columns) 
could  be: 


1 

2 

3 

4 

1 

2 

3 

4 

Mean 

4 

4 

+ 

+ 

Mean 

+ 

+ 

4 

4 

Ri 

- 

+ 

0 

0 

R. 

+ 

0 

0 

R2 

0 

0 

- 

+ 

or  R 

+ 

4 

-2 

0 

R3 

4 

+ 

- 

- 

R3 

+ 

4 

4 

-3 

The  interaction  of  row  and  column  contrasts  possibly  could  be  utilized  to  allocate 
the  symbols  in  the  latin  square. 

We  wish  to  illustrate  the  method  of  constructing  latin  squares  using 

orthogonal  polynomial  coefficients.  We  shall  first  consider  the  construction  of 

three  mutually  orthogonal  latin  squares  of  order  4  and  then  we  shall  consider  the 

construction  of  a  single  latin  square  of  order  6.  In  the  preceding  table  on 

orthogonal  polynomials  for  n  =  4  denote  all  combinations  with  a  plus  sign  as 

belonging  to  (R^l^I  and  those  with  a  minus  sign  as  belonging  to  (R^C^o  . 

Do  likewise  for  the  R,-C„  and  R  C  effects.  Then,  the  four  latin  square 

V  V  o  o 

symbols  are  obtained  as  follows: 

(rlCl)1,  (RQcQ)r  (R^^  =  0000  +  0101  +  1010  +  1111  =  a 

<rlclV  <rocq,o>(rccc)o=  000.  ♦  0100  +  1011  +  1110.  b 

<rlcDo’  «qcq>0,  (RcCc)j.  0010  +  0111  +  1000  +  1101  =  c 

(RL°L>0’  <Wl-  <rccc)0  =  0011  +  0110  +  1001  +  1100  =  D 
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This  results  in  the  following  latin  square  of  order  4 


A 

B 

C 

D 

B 

A 

D 

C 

C 

D 

A 

B 

D 

C 

B 

A 

Likewise,  if  we  use  the  following  polynomial  contrasts  we  obtain  the  two  mutually 
orthogonal  mates  of  the  above  square: 

(RLc  >r  <rqCc)0,  (rccl)0  =  0001  +  0110  +  1000  +  llll  =  a 

(RLCqV  ^rqccV  {Rccl\  ~  0010  +  0101  +  1011  +  1100  =  P 

(^Cglo,  <rqCc)0,  (r cCL)l  =  0011  +  0100  +  1010  +  1101  =  y 

(RLCq)0>  (RqCc)i,  <RCCL)0  =  0111  +  1001  +  1110  +  0000  =  6 

and 

(RLCC>1,  <RqClV  (Rccq,i  =  0011  +  0101  +  1000  +  1110  =  1 

(RLCc)r  (rqcl)0»  (RcCq)0  *  0001  +  0111  +  1010  +  1100  =  11 

(RLCc)0»  <RcVl  =  0000  +  0110  +  1011  +  1101  =  III 

(RLCc^o»  (RQClV  <rccq>0  s  0010  +  0100  +  1001  +  1111  s  IV 

The  above  results  in  the  following  two  latin  squares  of  order  4 


6 

a 

P 

Y 

V 

P 

a 

6 

a 

6 

Y 

P 

P 

6 

a 

~uT 

II 

IV 

I 

IV 

I 

III 

II 

I 

IV 

II 

III 

II 

III 

I 

IV 

The  above  method  of  constructing  mutually  orthogonal  latin  squares  using 
polynomial  coefficients  works  for  latin  squares  of  order  n  where  n  =  2P  .  We 
need  another  procedure  for  other  values  of  n  and  shall  now  construct  a  latin 
square  of  order  6  from  the  orthogonal  polynomial  coefficients  in  the  table  of 
single  degree  of  freedom  contrasts  for  36  combinations.  If  we  observe  only 


the  signs  of  contrasts  we  note  that  the  36  combinations  may  be  classified  into 
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six  sets  of  four  with  like  signs  and  two  additional  sets  of  six.  The  latter  two 
sets  will  be  used  to  build  up  the  six  sets  of  four  into  six  sets  of  six  as  follows 
where  all  combinations  with  a  plus  sign  go  in  the  one  level  and  all  those  with  a 
minus  sign  go  in  the  zero  level  (see  page  32): 


^R2C2V  *R3C3V  ^R4C4^0’  ^R5C5^0  +  2  from  ^RlGlV  ^R2C2  VR3C3  V  ^R4C4^1’  ^S^Vl 

<RaC2,0’  "Wo-  "Wi  "Wo  + 

"Wo-  ,r3c3)o>  "Wo- "Wo* 

^R2C2^0*  (R3C3,r  (R4C4,0’  ^5^1  +  2  fr°m  ^l^Vo*  (R2°2V  (R3C3)0>  (R4Cv4)l,(R5C5)0 
(R3C3V  (R4C4V  (R5C5}0  + 

<r2C2V  Wo’  (R4C4>o’  Wi* 

From  these  sets  we  obtain 


(12  +  21  +  34  +  43)  +  (00  +  55)  =  A 

(02  +  20+  35  +  53)  +  (11  +  44)  =  B 

(01  +  10  +45  +  54)  +  (22  +  33)  =  C 

(04  +  15  +40  +  51)  +  (23  +  32)  =  D 

(03  +  25  +30  +  52)  +  (14  +  41)  =  E 

(13  +  24  +31  +  42)  +  (05  +  50)  =  F 


This  results  in  the  following  latin  square  of  order  6: 


00 

A 

10 

C 

20 

B 

30 

E 

40 

D 

50 

F 

01 

C 

11 

B 

21 

A 

31 

F 

41 

E 

51 

D 

02 

B 

12 

A 

22 

C 

32 

D 

42 

F 

52 

E 

03 

E 

13 

F 

23 

D 

33 

C 

43 

A 

53 

B 

04 

D 

14 

E 

24 

F 

34 

A 

44 

B 

54 

C 

05 

F 

15 

D 

25 

E 

35 

B 

45 

C 

55 

A 

The  pair  of  treatments  in  the  second  set  of  parentheses,  e.  g.  (00  +  55),  was 
picked  from  the  set  of  six  in  such  a  manner  as  to  have  i  and  j  in  the  combina¬ 
tion  ij,  contain  0,  1,  2,  3,  4,  and  5  since  each  letter  must  appear  once  in 
each  r  -w  and  once  in  each  column. 

It  would  be  interesting  and  perhaps  enlightening  to  carry  out  the  above 
procedure  for  n  =  10  and  12  and  to  exhaustively  study  the  complete  set  of  35 
contrasts  for  n  =  6  . 
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V.  Group  Construction  of  Q(nt  t)  Sets 


V.  0.  Introduction 

The  construction  of  0(n,t)  sets  bastd  on  groups  and  their 

associated  mappings  such  as  automorphism,  complete  mapping,  and  orthomorphism 

is  the  oldest  and  still  the  most  popular  method  for  n  not  of  the  form  4t  +  2  . 

Euler  [1782]  implicitly  utilized  some  properties  of  finite  groups  of  order  2t  +  1 

and  4t  for  his  construction  of  0(2t+l,2)  and  0(4t,2)  sets,  respectively.  It 

was  MacNeish  [1922]  who,  for  the  first  time,  explicitly  (however,  not  rigorously) 

utilized  group  properties  for  his  construction  of  0(qm,  qm-l)  sets  and  0(n,  X.  ) 

fl  f2  fr 

sets,  where  q  is  a  prime,  m  is  a  positive  integer  and  if  n  =  q  q  ...  q 

*1  i2  *  r 

is  the  prime  power  decomposition  of  n  then  X.  =  min(q  ,  q  ,. . .  ,  q  )  -1  .  The 

l  w  r 

field  construction  of  0(qm,qm-l)  sets  found  independently  by  Bose  [1938] 
and  Stevens  [1939]  is  based  on  the  additive  group  of  GF(qm)  and  its  related 
cyclic  group  of  automorphisms.  The  0(n,n-l)  sets  for  n  =  3, 4, 5, 7, 8  and  9 
exhibited  by  Fisher  and  Yates  [1957]  are  based  on  cyclic  group  and  abelian  groups. 
Several  beautiful  applications  of  group  theory  to  the  existence  and  non-existence 
of  0(n,t)  sets  have  been  found  by  Mann  [1942,  1943,  1944]  .  The  0(12,5)  sets 
found  by  Johnson  ,et_  al.  [1961]  and  Bose  et  al.  [I960]  are  based  on  abelian 
groups  of  order  12.  Hedayat  [1969]  and  Hedayat  and  Federer  [1969]  have  found 
a  series  of  results  on  the  existence  and  non-existence  of  0(n,t)  sets  through  the 
group  theory  approach.  The  reader  interested  in  this  subject  will  find  the  fol¬ 
lowing  references  together  with  the  references  given  to  these  papers  very  useful: 
Paige  [1951],  Hall-Paige  [1955],  Singer  [I960],  Bruck  [  1951] ,  and  Sade  [1958]. 
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The  author  has  no  doubt  that  the  reader  can  find  many  more  interesting  papers 
directly  or  indirectly  related  to  this  rich  subject. 

V.  1.  Definitions  and  Notations 

There  are  several  forms  of  definitions  of  latin  squares  and  orthogonal 
latin  squares.  The  following  forms  are  useful  for  the  results  which  will  follow: 
Definition  V.  1. 1.  A  latin  square  of  order  n  on  an  n-set  ^  is  an  n  x  n  matrix 
whose  rows  and  columns  are  each  a  permutation  of  the  set  Every  latin  square 
of  order  n  may  therefore  be  identified  with  a  set  of  n  permutations  (p  ,p  ,. . .  ,p  ) 
where  pA  is  the  permutation  associated  with  the  ith  row. 

Definition  V.  1.  2  .  Let  Li  be  a  latin  square  of  order  n  on  an  n-set  ^  , 
i  =  1, 2, . . .  ,t  .  Then,  the  set  S  =  {Lj,  L2,  . . . ,  L  }  is  said  to  be  a  mutually 
orthogonal  set  of  t  latin  squares  if  the  projection  of  the  superimposed  form  of 
the  t  latin  squares  on  any  two  n-sets  ^  and  Zj>  i*J>  forms  a  permutation 
of  the  cartesian  product  set  of  ^  and  Zj  •  Such  a  set  is  denoted  as  an 
0(n,t)  set.  (See  also  definitions  I.  2  and  1. 3 . ) 

Definition  V.  1.  3 .  If  Lj  =  (P^,  P12,  . . . ,  PJn)  and  L2  =  (P21,  P22>  . . .  ,  P2n) 
are  two  latin  squares  of  order  n  on  an  n-set  Yj>  then  we  may  define  L^L2  to 

be  L3  =  (p11P2i>  P12P22*  *  *  *  *  PlnP2n^soe  de^inition  !)•  11116  generalization 
to  the  product  of  t  >  2  latin  squares  follows  immediately. 

V.  2  .  Construction  of  O(n.t)  Sets  Based  on  a  Group 

We  shall  divide  the  problem  into  three  parts  based  on  whether  n  is  a 
prime,  or  a  mixture  of  prime  powers.  The  proof  of  the  subsequent  results  can  be 
found  in  the  references  related  to  this  section. 
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V.  2. 1.  n  =  q  a  prime .  Recall  that  any  prime  ordered  group  is  cyclic. 

Theorem  V.  2. 1. 1.  Let  G  =  {P. ,  P- , . . .  ,  P  }  be  a  cyclic  permutation  group  of  de- 

1  M  CJ 

gree  q  and  order  q  .  Then.  =  {Lj,  L2 . L  j}  is  an  0(q,  q-1)  set. 

where  =  (pj,  P*,  ... ,  P*)  . 

Demonstration  V.  2. 1, 1.  Let  q  =  5  .  Select  any  arbitrary  generator  such  as 
^35  2  14^  which  generates  a  cyclic  permutation  group  G  and,  hence,  a  latin 
square  L^  .  Then, 


3 

5 

2 

1 

4 

2 

4 

5 

3 

1 

5 

1 

4 

2 

3 

4 

3 

1 

5 

2 

1 

2 

3 

4 

5 

•  L2 


2 

4 

5 

3 

T 

4 

3 

1 

5 

2 

3 

5 

2 

1 

4 

5 

1 

4 

2 

3 

1 

2 

3 

4 

T 

-  L3 


5 

1 

4 

T 

T 

3 

5 

2 

i 

4 

4 

3 

1 

5 

2 

2 

4 

5 

3 

1 

1 

2 

3 

4 

T 

-  V 


4 

3 

1 

5 

~T 

5 

1 

4 

2 

3 

2 

4 

5 

3 

1 

3 

5 

2 

1 

4 

1 

2 

3 

4 

5 

For  those  who  do  not  like  to  work  with  permutation  groups  we  present  the  following 
theorem: 

Theorem  V.  2. 1.  2.  igl  L(r)  be  an  n  X  n  square  with  ri  +  j  (mod  q)  in  its  (i,j)th 
cell,  i,  j  =  0,1,...  ,q-l  .  Then.  S12  =  {L(l),  L(2).  . . .  ,  L(q-l)}  is  an  0(q,  q-1) 
set  if  q  is  a  prime. 

Demonstration  V.  2. 1.  2.  Let  q  =  5;  then, 


L(l)  = 


0 

1 

T 

T 

4 

1 

2 

3 

4 

0 

2 

3 

4 

0 

1 

3 

4 

0 

1 

2 

4 

0 

1 

2 

3 

,  L(2)  = 


o" 

1 

2 

3 

4 

2 

3 

4 

0 

1 

4 

0 

1 

2 

3 

1 

2 

3 

4 

0 

3 

4 

0 

1 

2 

,  L(3)  = 


0 

1 

T 

3 

4 

3 

4 

0 

1 

2 

1 

2 

3 

4 

0 

4 

0 

1 

2 

3 

2 

3 

4 

0 

1 

,L(4)  = 


0 

1 

2 

3 

4 

4 

0 

1 

2 

3 

3 

4 

0 

1 

2 

2 

3 

4 

0 

1 

1 

2 

3 

4 

0 
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Note  that  L(  1)  in  theorem  V.  2. 1.  2  is  based  on  the  cyclic  permutation 

gioup  generated  by  2  3  * '  *  o  ^  and  L(i)  =  1  =  2>3 . q_1  *  Hence 

theorem  V.  2. 1.  2  is  a  special  case  of  theorem  V.  2. 1. 1  . 

V.  2.  2.  n  =  qm  where  q  is  a  prime  and  m  any  positive  integer.  Note  that 

this  case  in  particular  for  m  =  1  includes  case  1.  We  shall  present  three 

theorems  for  this  case.  The  first  two  are  based  on  cyclic  groups  and  the  third 

one  is  based  on  any  group  which  admits  an  automorphism  of  order  t  . 

Theorem  V.  2.  2. 1.  Let  G  =  {P. ,  P? ,  . . . ,  P  }  be  a  cyclic  permutation  group  of 

degree  n  and  order  n  .  Then.  S21  =  {Lp L2,. . .  ,  L  }  is  an  0(n, \  )  set 

where  n  =  qm  ,  \  =  q  - 1,  and  =  (p{ ,  P*,  •  •  • ,  P* )  • 

2 

Demonstration  V.  2.  2. 1.  Let  n  =  3  =  9  .  Select  any  arbitrary  generator  such 

as  (i  !  t  1  i  ?  !  !  !)  which  generates  a  cyclic  permutation  group  G  and  hence, 

1  O  107c, 

a  latin  square  L  .  Then,  since  X.  =  2, 


and  L2  = 


is  an  0(9 , 2)  set. 

Conjecture.  The  set  is  orthogonally  locked,  meaning  that  there  does  not 

-!<  ..  * 

exist  a  latin  square  L  such  that  S2j  U  {L  }  is  an  0(n,\  +1)  set  . 


5 

1 

6 

3 

7 

~ 

9 

2 

~ 

7 

5 

8 

6 

9 

2 

4 

1 

3 

9 

7 

2 

8 

4 

1 

3 

5 

6 

4 

9 

1 

2 

3 

5 

6 

7 

8 

3 

4 

5 

1 

6 

7 

8 

9 

2 

6 

3 

7 

5 

8 

9 

2 

4 

1 

8 

6 

9 

7 

2 

4 

1 

3 

5 

2 

8 

4 

9 

1 

3 

5 

6 

7 

1 

2 

3 

4 

5 

6 

7 

8 

9 

3 

4 

5 

1 

6 

7 

8 

9 

2 

5 

1 

6 

3 

7 

8 

9 

2 

4 

6 

3 

7 

5 

8 

9 

2 

4 

1 

7 

5 

8 

6 

9 

2 

4 

1 

3 

8 

6 

9 

7 

2 

4 

1 

3 

5 

9 

7 

2 

8 

4 

1 

3 

5 

6 

2 

8 

4 

9 

1 

3 

5 

6 

7 

4 

9 

1 

2 

3 

5 

6 

7 

8 

1 

2 

3 

4 

5 

6 

7 

8 

9 
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Note  that  for  n  even  this  conjecture  is  correct  since  any  latin  square 
of  even  order  based  on  cyclic  permutation  group  is  orthogonally  mateless. 

An  analogous  theorem  to  theorem  V.  2. 1.  2  for  this  case  is: 

Theorem  V.  2.  2.  2.  Let  L(r)  be  an  n  X  n  square  with  ri  +  j  (mod  n)  in  its 
(i,  j)  cell,  i  =  0,1,2,...  ,n-l  .  Then  =  {L(l ) ,  L(2),  .  .  .  ,  L(  X  )}  is  an 
0(n,\)  set  if  n  =  qm  and  X  =  q-1  . 

2 

Demonstration  V.  2.  2.  2.  Let  n  =  q  =  3  then, 


L(l)  = 


0 

1 

2 

3 

4 

5 

6 

7 

8 

1 

2 

3 

4 

5 

6 

7 

8 

0 

2 

3 

4 

5 

6 

7 

8 

3 

4 

5 

6 

D 

8 

0 

1 

2 

4 

5 

6 

7 

8 

0 

D 

3 

5 

6 

7 

8 

0 

1 

2 

3 

4 

6 

7 

8 

0 

1 

2 

3 

4 

5 

7 

8 

0 

1 

2 

3 

4 

5 

6 

8 

0 

1 

2 

3 

4 

5 

6 

7 

and  L(2) 


0 

1 

2 

3 

4 

5 

6 

7 

8 

2 

3 

. 

4 

5 

6 

7 

8 

0 

1 

4 

5 

6 

7 

8 

0 

1 

2 

3 

6 

7 

8 

0 

1 

2 

3 

4 

5 

8 

0 

1 

2 

3 

4 

5 

6 

7 

1 

2 

3 

4 

5 

6 

7 

8 

0 

3 

4 

5 

6 

7 

8 

0 

1 

2 

5 

6 

7 

8 

0 

1 

2 

3 

4 

7 

8 

0 

1 

2 

3 

4 

5 

6 

is  an  0(9,2)  set.  Note  that  theorem  V.  2.  2.  2  is  a  special  case  of  theorem  V.  2 . 2 . 1 
viz.  ,  L(l)  is  based  on  the  cyclic  permutation  group  generated  by  (^  ^  n  ) 

and  L(i)  =  L1  ( 1 ) ,  i  =  2, . . .  ,  X  . 


Theorem  V.  2.  2.  3.  Let  G  =  {a^  =  e  the  identity,  a^,  . . . 
order  n  and  a  an  automorphism  on  G  such  that  a^a^)  + 
aj  *  e  . 

1)  S  =  {Lj,  L^,  .  .  ,  Lt)  is  an  0(n,  t)  s^i,  where 


,  a^}  be  a  group  of 
a.,  1  <  i  <  t  . 
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e 

a2 

•  •  • 

an 

«‘(a2l 

«\a2)a2 

•  •  • 

or^aja 

2  n 

^(a3) 

*‘<a3)a2 

•  •  • 

a4(a  )a 

3  n 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

al( a  ) 

^(a  )a 

•  •  • 

aA(a  )a 

n 

n  2 

n  n 

i  l,2,...,t  • 

2)  If  in  particular  t  =  n-1,  then  one  can  simplify  the  construction  of  an 
0(n,  n-1)  set  from  the  following  key  latin  square  by  a  cyclic  permutation  of  its 
last  n-1  rows. 


e 

a(x) 

2.  . 

a  (x) 

•  •  • 

t,  V 
a  (X) 

«(x) 

a(x)a(x) 

a(x)a2(x) 

•  •  • 

a(x)at(x) 

2.  » 
a  (x) 

2 

a  (x)a(x) 

2  2 
a  (x)a  (x) 

•  •  • 

ar^(x)ar*(x) 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

^(x) 

a  (x)a(x) 

t.  2 
a  (x)a  (x) 

•  •  • 

at(x)at(x) 

for  anv  x  _i&  G  except  the  identity  element. 

We  see,  therefore,  that  by  means  of  theorem  V.  2.  2.  3  one  can  construct 
an  0(n,t)  set  if  we  can  find  a  group  G  and  an  automorphism  a  of  order  t  . 

In  particular,  if  t  =  n-1  the  whole  task  of  construction  reduces  to  the  construc¬ 
tion  of  Lq  as  described  above.  If  n  =  qm  then  because  every  elementary 
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abelian  q-group  G  of  order  n  admits  an  automorphism  a  of  order  n-1,  we 
can  construct  an  0(qm,qm-l)  set  based  on  G  and  a  .  Here  we  present  a 
general  method  of  constructing  such  an  automorphism  for  any  n  =  qm  .  In  par¬ 
ticular,  we  shall  exhibit  such  an  automorphism  for  the  following  n  : 

to  n 

n  =  2  ,  m  =  2 ,  3 ,  . . .  ,  9 

n  =  3m,  m  =  2,  3,  . . .  ,  6 

n  =  5m,  m  =  2,  3,  4 

n  =  7  ,  m  =  2,  3 

n  =  ll2,  132,  172,  192,  232,  292,  and  312  . 

This  will  then  perhaps  be  the  largest  table  that  has  ever  been  produced  so  far 
for  0(n,n-l)  sets. 

Note  that  there  is  no  loss  of  generality  if  we  limit  ourselves  to  the  follow¬ 
ing  elementary  abelian  q-group  of  order  n  =  qm  . 

?!< 

G  —  { (bj  b^  •  •  •  b^),  bj  -  0,1,  2,  ,  •  •  ,  q— lj  j  -  1,2,...,  in}  . 

>;< 

The  binary  operation  on  G  is  addition  mod  q  componentwise,  viz. .  (b^  b2  . . . 

b  )  +  (bl  b'  . . .  b'  )  =  (c.  c_  . . .  c  )  where  c.  =  b.  +  b!  (mod  q)  .  Note  that 
m  1  2  m  12  m  111 

>!< 

the  elements  of  G  are  simply  the  treatment  combinations  of  m  factors  each 
at  q  levels.  The  reason  why  we  have  chosen  this  particular  elementary  abelian 
q-group  is  that  it  has  a  well-known  structure  to  those  who  are  concerned  with  ex- 
periment  design  construction.  Note  also  that  G  is  the  direct  product  of  m 
Galois  fields,  each  of  order  q  . 

The  generator  set  for  every  elementary  abelian  q-group  of  order  qm  con¬ 
sists  of  m  elements,  and  for  uniformity,  we  may  choose  the  following  ordered 
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generator  set  for  G  . 

g  =  {(100  ...  0),  (01,00  ...  0),  ...  (00  ...,  010),  (00  ...  01)}  . 

>:< 

Note  that  the  structure  of  every  automorphism  a  on  G  is  completely 

5j; 

defined  if  we  know  the  image  of  each  element  of  g  under  a  .  G  is  a  vector 

space  of  dimension  m  over  GF(q). 

Before  proceeding  further  we  need  the  following  well-known  result: 

Theorem  V.  2.  2.  4.  Let  G  be  an  elementary  abelian  q-qroup  of  order  n  =  qm  . 

Then,  Automorphism  group  of  G  is  isomorphic  to  the  (multiplicative)  group  of  all 

non-singular  m  X  m  matrices  with  entries  in  the  field  of  integers  mod  q  . 

The  construction  of  an  automorphism  of  order  n-1  for  G  is  equivalent 

1  ^ 

to  the  construction  of  an  m  X  m  matrix  A  such  that  A  =  I  but  A  *  I  if  t 

is  not  a  multiple  of  n-1  over  the  field  of  integers  mod  q  . 

We  know  from  linear  algebra  that  if  4»  is  a  linear  map  on  a  vector  space 

V  and  if  x  e  V  such  that  x  *  0  but  <{>(x)  =  x,  then  1  is  an  eigenvalue  of  <j>. 

f  s  s 

Moreover,  if  \  ,  . . .  ,  X.  }  is  the  set  of  eigenvalues  of  <f>,  then 

1  4.  t  1  '£ 

S  S 

. . .  ,  \  }  is  the  set  of  eigenvalues  of  «J>  .  Therefore,  for  our  problem  we  must 

find  a  linear  map  on  G  with  a  set  of  eigenvalues  \  having  the  property  that 

s  n-1 

for  each  i,  X.  *  1  (mod  q )  for  all  s  =  1, 2, . . .  ,n-2  and  X. ,  =  1  .  To  do  so 

let  F  be  a  GF(qm)  and  let  (3  be  a  generator  of  the  multiplicative  cyclic  group 

of  GF(qm),  i.e.  (3**1,  i  =  1, 2, . . .  ,n-2  while  (3n  1  =  1  .  Let  f(x)  be  a 

monic  irreducible  polynomial  over  GF(q)  for  (3  .  Note  that  f(x)  has  degree  m  . 

(3  is  sometimes  called  a  primitive  root  or  mark  of  F  .  Now,  if  we  let  A  be  the 

companion  matrix  for  (3,  then  it  is  easy  to  see  that  A  has  the  desired  property. 
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Example 


Let  us  find  an  automorphism  of  order  3  for  G  =  {(00),  (01),  (10),  (11)}. 
It  is  sufficient,  by  previous  arguments ,  to  find  a  2X2  matrix  A  of  order  3 
over  the  field  of  integer  mod  2  .  Let  GF(22)  =  (o,  1,  f3,  (3+  1}  with  the  follow¬ 
ing  addition  (+)  and  multiplication  (•  )  tables 


he 

1 

p 

p+ 1 

0 

0 

1 

p 

p+ 1 

1 

1 

0 

p+ 1 

p 

p 

p 

p  + 1 

0 

1 

p  + 1 

p  +  1 

p 

1 

0 

• 

0 

1 

p 

p+ 1 

0 

0 

0 

0 

0 

1 

0 

1 

p 

p+  1 

p 

0 

p 

P-+  1 

1 

p  + 1 

0 

p  +  1 

1 

p 

Note  that  (3  is  a  primitive  root  for  GF(22)  and  f(x)  =  x2  +  x  +  1  is  a  monic 
irreducible  polynomial  for  (3,  since  f((3)  =  0(mod2)  .  The  companion  matrix 
associated  with  f(x)  is 


A  = 


0  1 
1  1 


As  a  check 


1 

1  1 

3 

1  2 

■»  mm 

1  0 

= 

over  GF(2),  A  = 

= 

2_ 

_i  0_ 

_0  1_ 

-0  L 

over  GF(2)  . 

Let  us  now  determine  the  image  of  the  ordered  generator  set  g  =  {(10),  (01)} 
under  A  . 
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o  (10)  t  i  (01) 


Ag  = 


o  r 

"(10)' 

= 

i  i 

(01)_ 

1  (10)  4  1  (01) 


"(01)' 

= 

Therefore,  A(10)  =  (01),  A(01)  =  (11),  and  since  (11)  «  (10)  +  (01),  (00)  =  2(10) 

+  2(01),  we  have  A(ll)  =  (10),  A(00)  =  (00)  . 

Now,  we  have  a  group  G  of  order  4  and  an  automorphism  of  order  3 
on  G  .  We  can  now  construct  an  0(4,3)  set  .  Since  e  =  (00),  and  If  we  let 
x  =  (10)  In  theorem  V.  2.  2.  3,  we  obtain: 


(00) 

A(10) 

AZ(10) 

a3(10) 

A(10) 

A(10)  A(10) 

A(10)  A2 (10) 

A(10)  A3(10) 

A2  (10) 

A2(10)A(10) 

A2(10)A2(10) 

A2(10)A3(10) 

A3(10) 

A3(10)A(10) 

A3(10)A2(10) 

A3(10)A3(10) 

(00) 

(01) 

(11) 

(10) 

(01) 

(00) 

(10) 

(11) 

(11) 

(10) 

(00) 

(01) 

(10) 

(11) 

(01) 

(00) 

The  other  two  la  tin  squares  are  obtained  by  a  cyclic  permutation  of  the  last  three 
rows  of  Lq  .  Thus, 


(00) 

(01) 

(11) 

(10) 

(11) 

(10) 

(00) 

(01) 

(10) 

(11) 

(01) 

(00) 

(01) 

(00) 

(10) 

(ii) 

(00) 

(01) 

(11) 

(10) 

(10) 

(11) 

(01) 

(00) 

(01) 

(00) 

(10) 

(11) 

(ID 

(10) 

(00) 

(01) 
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To  simplify  the  notation  we  set  (00)  =  1,  (01)  =  2,  (11)  =  3,  (10)  =  4  to  obtain: 


1 

2 

3 

4 

2 

1 

4 

3 

3 

4 

1 

2 

4 

3 

2 

1 

1 

2 

3 

4 

4 

3 

2 

1 

2 

1 

4 

3 

3 

4 

1 

2 

and  L 


1 

2 

3 

4 

3 

4 

1 

2 

4 

3 

2 

1 

2 

1 

4 

3 

We  are  now  ready  to  exhibit  a  generating  matrix  of  order  n-1  =  q111 
with  entries  from  GF(q)  for  those  n  promised  before. 


n 


Generator 


Order 


n 


Generator 


Order 


l” 

0 

1 

0 

1 

X 

1 

3 

23 

0 

0 

1 

— 

— 

1 

0 

1 

— 

— 

— 

— * 

— 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

5 

0 

0 

1 

0 

0 

0 

0 

0 

1 

15 

2 

0 

0 

0 

1 

0 

1 

1 

0 

0 

mm 

0 

0 

0 

0 

1 

1 

0 

1 

0 

0 

IB 

f— 

-1 

_ 

— 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

63 

7 

2 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

1 

1 

0 

0 

0 

0^ 

1 

1 

0 

0 

0 

0 

0 

L 

J 

r 

— 

0  1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0  0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

255 

29 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

1 

1 

1 

0 

0 

0  _ 

0 

0 

□ 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

0 

0 

0 

1 

— 

7 


31 


127 


511 
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n 


3 


7 


Generator 


Order 


n 


Gen  erator 


0  1  0 
0  0  1 
1  1  0 


"o 

1 

0 

0 

"o 

1 

0 

0 

0 

0 

0 

1 

0 

80 

,5 

0 

0 

1 

0 

0 

3 

1 

0 

0 

0 

0 

1 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

1 

1 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

r\ 

\J 

1 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

i 

1 

1 

0 

0 

0 

0 

0  10  0 
0  0  10 
0  0  0  1 
3  3  0  3 


728 


624 


0  1  0 
0  0  1 
2  2  0 


342 


288 

528 

840 


0  1  0 
0  0  1 
2  2  0 


To  shed  more  light  on  the  given  procedure  we  go  through  another  example. 
3 

n  =  2  .  Then 


Order 

26 

242 

124 

48 

120 

168 

360 

728 

960 

Let 
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g"  =  {(000),  (001),  (010),  (Oil),  (100),  (101),  (110),  (111)} 

"o  1  0 

g  =  {(100),  (010),  (001)}  and  A  =  0  0  1 

1  0  1_ 

"o  1  o’]  p (100)1  r (010) 

Ag  =  001  (010)  =  (001) 

1  0  1  (001 )_  (101) 

Let  x  in  theorem  V.  2.  2.  3  be  (100).  Then, 

A(100)  =  (010)  , 

A2(100)=  (001)  , 

3 

A  (100)=  (101)  , 

A4(100)=  (111)  , 

A5(100)=  (110)  , 

A6(100)=  (Oil)  ,  and 
A7(100)=  (100)  . 

Therefore,  we  obtain  LQ  as  follows: 


(000) 

(010) 

(001) 

(101) 

(111) 

(110) 

(Oil) 

(100) 

(010) 

(000) 

(Oil) 

(111) 

(101) 

(100) 

(001) 

(110) 

(001) 

(011) 

(000) 

(100) 

(110) 

(111) 

(010) 

(101) 

(101) 

(111) 

(100) 

(000) 

(010) 

(Oil) 

(110) 

(001) 

(111) 

(101) 

(110) 

(010) 

(000) 

(001) 

(100) 

(Oil) 

(110) 

(100) 

(111) 

(Oil) 

(001) 

(000) 

(101) 

(010) 

(Oil) 

(001) 

(010) 

(110) 

(100) 

(101) 

(000) 

(111) 

(100) 

(110) 

(101) 

(Oil) 

(Oil) 

(010) 

(111) 

(000) 

Setting  (000)  =  1,  (010)  =  2,  (001)  =  3,  (101)  =  4,  (111)  =  5,  (110)  =  6  , 

(0  11)  =  7,  (100)  =  8,  then  LQ  in  a  compact  form  will  be: 
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Now,  we  can  derive  Lj,  L_»  •  •  •  >  from  LQ  by  a  cyclic  permutation  of  the 
last  7  rows  of  LQ,  for  example, 


and  so  on.  Note  the  way  is  derived  from  LQ  :  except  for  the  first  row  of 

Lq  and  Lj,  which  are  identical,  the  i^  row  of  LQ  becomes  the  (i  +  1)^  row 

of  L^,  and  the  last  row  of  LQ  becomes  the  second  row  of  .  In  general 

is  derived  from  L.  ,  in  the  same  fashion  as  L,  is  derived  from  Ln  . 

j-1  I  0 


m  m^  m^ 

V.  2.  3.  n  =  q1  1  q2  . . .  qr  r,  where  qA  is  a  prime  such  that  q,  *  q,  if  i  *  j 


and  m,  is  a  positive  integer,  i  =  1,2,...  ,r  . 
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m.  m_  m 

12  r 

Theorem  V.  2.  3,1.  Let  n  =  q^  ...  qr  be  the  prime  power  decomposition 

ml 

of  n  .  Then,  there  exists  an  0(n,,y)  set  based  on  a  group,  where  y  =  min(q^  , 

m  m 

2  1 

^2  >  •  •  •  >  f  • 

m. 


Construction.  Let  n^  =  q,  .  Then,  by  the  method  of  theorem  V.  2.  3  construct 

®n  0(n  ,  nj-1)  set  S.  =  {L^,  Li2,  . . .  ,  L.n  _]) ,  i=l,2,...,r  .  Now,  let 

Sj  ~  (L^,  L^ >  •  •  •  >  L^} ,  i  =  1,  2,  . . .  ,  r  .  Then,  H  -  {A^,  ,  •  •  •  ,  A  }  is 

an  0(n,  y)  set  where  A,  =  L](®L  ®...  $L  ,  .  ®  denotes  the  Kronecker  product 

J  c]  rj 

operation. 

2 

Demonstration  V.  2.  3. 1 .  Let  n  =  12  =  2  •  3  .  Then,  -y  =  2, 


V  Lu  = 


1 

~T~ 

3 

2 

3 

1 

3 

1 

2 

-  Lu  = 


1 

2 

3 

3 

1 

2 

2 

3 

1 

s 


2 


L 


1 

2 

~3 

T 

2 

1 

4 

3 

3 

4 

1 

2 

4 

3 

2 

1 

T 

2 

T 

4 

4 

3 

2 

1 

2 

1 

4 

3 

3 

4 

1 

2 

[Tj 

"2" 

T 

4 

3 

4 

1 

2 

4 

3 

2 

1 

2 

1 

4 

3 

S1  *tLll'L12}  and  S2  SM  (L21>L22} 
verify  that 


Then,  the  reader  can  easily 


H  =  {Aj  =  L^  $  L21,  A2  =  L12®  L22  } 

is  an  0(12,  2)  set. 

Remark.  Let  n  and  y  be  the  same  as  in  theorem  V.  2.  3.1.  Then  it  can  be  shown 
that  automorphism  method  fails  to  produce  more  than  y  mutually  orthogonal  latin 


#1030 


-47- 


squares.  We  shortly  show  that  this  inherent  defect  is  due  to  the  mapping  not 
to  the  group  structure. 

Definition  V.  2. 1.  Consider  for  each  positive  integer  n  an  abstract  group  G  of 
order  n  with  binary  operation  *  .  Let  D  be  the  collection  of  all  one-to-one 
mappings  of  G  into  itself.  Then  two  maps  a  and  in  D  are  said  to  be 
orthogonal  if  for  any  g  in  G  . 

(a-  Z)  *  (4j  Z)  1  =  g 

has  a  unique  solution  Z  in  G  .  In  particular  if  <r  is  an  identity  map  then 
is  said  to  be  an  orthomorphism  map.  A  t-subset  of  11  is  said  to  be  a  mutually 
orthogonal  set  if  every  two  maps  in  this  t-subset  are  orthogonal. 

Let  L(* )  be  an  n  X  n  square.  We  make  a  one-to-one  correspondence 
between  the  rows  of  L(*  )  and  the  elements  of  G  .  Thus,  by  row  x  we  shall 
mean  the  row  corresponding  to  the  element  x  in  G  .  Similarly  we  make  a  one- 
to-one  correspondence  between  the  columns  of  L(*  )  and  the  elements  of  G  . 

The  cell  of  L(- )  which  occurs  in  the  intersection  of  row  x  and  column  y  is 
called  the  cell  (x,y)  . 

Theorem  V.  2.  3.  2.  Let  <j  be  in  11  .  Put  in  the  cell  (x,y)  of  L(* )  the  element 
(o-x)  *  y  of  G  .  Call  the  resulting  square  L(<r  )  .  Then  L(o- )  is  a  latin  square 
of  order  n  on  G  .  Moreover  if  {o-  ^ ,  <r  . . .,  o- is  a  set  of  t  mutually  ortho¬ 
gonal  maps  then  {L(o-j),  . . .  ,  L(<r  )}  is  an  C(n,t)  set. 

Demonstration  V.  2.  3,  2.  Let  G  =  {0,  1,  2}  with  the  binary  operation  x  +  x  = 

1  M 

x3(mod  3),  xi  in  G  .  Then  the  maps  <r  and  with  the  following  definitions 
are  orthogonal. 
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cr  (0)  =  0 

<7(1)  =  1 
«r(2)  =  2 


<M0)  =  o 
4>(D  =  2 

4j(2)  =  1  . 


The  corresponding  latin  squares  to  o-  and  ijj  are: 


0 

1 

2 

1 

2 

0 

2 

0 

1 

L(4>)  = 


0 

1 

2 

2 

0 

1 

1 

2 

0 

which  are  orthogonal. 


V.  3.  Construction  of  0(n,t)  sets  based  on  t  different  groups  of  order  n 

Up  to  now  we  have  been  concerned  with  the  construction  of  0(n,t)  sets 
using  a  group  of  order  n  which  admits  certain  mappings.  In  this  section  we  want 
to  show  that  for  some  n's  and  t's  one  can  construct  0(n,t)  sets  based  on  t 
different  groups  each  of  order  n  .  This  approach  proved  useful  because  it  lead 
to  the  construction  of  an  0(15,3)  set.  We  should  mention  that  our  motivation 
to  search  along  these  lines  has  stemmed  from  the  following  theorem,  with  a  nega¬ 
tive  flavor,  proved  by  Mann  [1944]  . 

Theorem  V.  3. 1.  It  is  impossible  to  construct  an  0(5,2)  set  based  on  two  dif¬ 
ferent  permutation  groups. 

For  a  while  we  thought  that  this  theorem  might  be  true  for  all  other  orders. 
However,  it  was  found  that,  fortunately,  this  is  not  the  case  as  the  following  two 
theorems  show: 

Theorem  V.  3.  2.  It  is  possible  to  construct  0(7,2)  sets  based  on  two  different 
cyclic  permutation  groups  of  order  7  . 
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Proofs  (Lj,  L2)  is  an  0(7,2)  set  where 


1 

2 

3 

4 

5 

6 

7 

3 

7 

6 

1 

4 

2 

5 

6 

5 

2 

3 

1 

7 

4 

2 

4 

7 

6 

3 

5 

1 

7 

1 

5 

2 

6 

4 

3 

5 

3 

4 

7 

2 

1 

6 

4 

6 

1 

5 

7 

3 

2 

T 

2 

3 

~ 

5 

6 

7 

2 

3 

4 

5 

6 

7 

1 

3 

4 

5 

6 

7 

1 

2 

4 

5 

6 

7 

1 

2 

3 

5 

6 

7 

1 

2 

3 

4 

6 

7 

1 

2 

3 

4 

5 

7 

1 

2 

3 

4 

5 

6 

Lj  and  L  are  based  on  two  different  permutation  groups  as  can  easily  be  seen 

from  the  different  structure  of  their  rows.  To  be  specific  is  based  on  the 

1  Z  3  4  5  6  7 

cyclic  permutation  group  generated  by  (3  7  6  j  4  2  5^  and  L2  is  based  on  the 
cyclic  permutation  group  generated  by  (3  3  4  5  6  7  1>  *  Note  t*iat>  since  Li  and 
L2  are  based  on  cyclic  permutation  groups,  then  by  theorem  V.  2.1. 1  {L^}  and 
{L2}  can  be  embedded  in  0(7,6)  sets.  However,  whether  or  not  {L^,  L2)  can 
be  embedded  in  a  larger  set  is  an  open  problem. 

Theorem  V.  3.  3.  It  is  possible  to  construct  0(15,  3)  sets  based  on  three  dif¬ 
ferent  cyclic  permutation  groups  of  order  15. 

We  remind  the  reader  that  every  group  of  order  15  is  cyclic. 

Proof:  {Lj,  L2,  is  an  0(15,  3)  set  where 
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0 

i 

2 

3 

4 

~T 

3 

4 

5 

6 

7 

U 

5 

6 

7 

8 

9 

Ei 

■a 

8 

9 

10 

1 1 

8 

9 

10 

1 1 

12 

13 

m 

11 

12 

13 

14 

0 

13 

14 

0 

1 

2 

El 

1 

2 

3 

4 

u 

2 

3 

■a 

5 

D 

MM 

5 

6 

mm 

8 

B 

6 

7 

8 

9 

10 

§1 

8 

9 

10 

1 1 

12 

B 

10 

11 

12 

13 

14 

ID 

12 

13 

14 

0 

1 

14 

0 

1 

2 

generated  by 


fo  1  2  3  4  5 

\2  3  4  5  6  7 


H 

2 

3 

4 

5 

rr 

11 

10 

MEM 

9 

14 

2 

4 

6 

3 

12 

■a 

10 

9 

13 

mm 

8 

MM 

3 

5 

6 

ni 

S 

9 

7 

14 

13 

i 

9 

3 

6 

12 

5 

1 1 

u 

u 

13 

8 

14 

2 

B 

mm 

5 

0 

12 

10 

L<l 

13 

14 

1 

8 

4 

5 

12 

1 1 

0 

9 

m 

14 

8 

2 

1 

3 

12 

0 

10 

1  1 

7 

ID 

8 

1 

4 

2 

6 

LL 

0 

1  1 

9 

10 

13 

generated  by 


0  1  2  3  4  5 
1  11  10  7  9  14 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

7 

4 

8 

9 

1  1 

2 

5 

12 

13 

10 

0 

14 

1 

3 

12 

1  1 

13 

10 
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G 


1  2  3  4  5  6  7  8  9  10  11  12  13 

4  8  9  11  2  5  12  13  10  0  14  1  3 


li)  • 


Whether  or  not  {l^,  L 2,  L^}  can  be  embedded  in  an  0(15,  t ),  t  >3,  set 
is  an  open  problem  (see  Hedayat  f  1970 ] ). 

V.  4.  Concluding  Remark 

Johnson  et  al.  [1961]  and  Bose  etal.  [I960]  independently  found,  by  an 
electronic  computer,  five  mutually  orthogonal  latin  squares  by  first  finding  five 
mutually  orthogonal  maps  for  an  abelian  group  of  order  12.  The  0(12,  5)  set 
exhibited  below  is  the  set  found  by  Johnson  et  al.  [1961]  .  Note  that  the  top 
square  is  obtained,  after  a  proper  renaming,  as  the  direct  product  of  a  latin  square 
of  order  2  and  a  cyclic  latin  square  of  order  6  being  both  orthogonally  mateless. 
Moreover,  every  other  square  is  obtained  by  proper  row  permutations,  determined 
by  an  orthomorphism,  from  the  top  square. 

Final  Remark.  The  group  method  fails  to  produce  an  0(n,t)  set,  t  >  2  for  any 

n  of  the  form  4t  +  2  .  This  is  so  because  the  Cayley  table  of  any  group  of  order 
n  =  4t  +  2,  which  is  a  latin  square  of  order  n,  is  orthogonally  mateless. 
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VI.  Projecting  Diagonals  Construction  of  Q(n,  t)  Sets 


A  very  simple  procedure  (sort  of  the  "man-on- the-street"  approach)  of 
1 

constructing  balanced  incomplete  block  and  partially  balanced  incomplete  block 
2 

designs  for  v  =  k  items  in  incomplete  blocks  of  size  k  has  been  utilized 
since  the  late  1940' s  by  the  author  and  has  its  counterpart  in  constructing  0(n,t) 
sets.  First  we  shall  illustrate  its  use  in  incomplete  block  experiment  design 
construction,  and  then  we  show  how  it  applies  to  the  construction  on  0(n,t)  set. 
The  theoretical  basis  for  this  method  may  be  derived  directly  from  the  preceding 
section. 

The  procedure  becomes  apparent  through  an  example.  Suppose  that  v  =  9 
and  k  =  3  .  After  writing  the  first  square  as  illustrated  below,  take  successive 
diagonals  of  the  preceding  square  and  use  them  to  form  the  incomplete  blocks  of 
a  square,  thus: 


Square  1  Square  2  Square  3 


Square  4 


1 

4 

7 

2 

5 

8 

3 

6 

9 

As  we  have  noted  this  is  a  resolvable  balanced  incomplete  design  with  the  para- 
2 

meters  v  =  9  =  k  ,k=3,  r  =  4  =  k  +  1,  b  =  12  =  k(k+l),  and  X.  =  1,  where 
the  rows  of  the  above  squares  form  the  incomplete  blocks. 

2 

To  form  a  partially  balanced  incomplete  block  design  for  v  =  k  in  in¬ 
complete  blocks  of  size  k  one  may  use  any  2,  any  3,  ...  ,  any  k  arrangements 
(or  squares).  To  illustrate  the  formation  of  a  partially  balanced  incomplete  block 
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design  for  v  =  6  =  k(k-l),  r  =  2,3,  . .  .,k,  and  k'  =  k-1  =  2  simply  delete  the 

last  set  of  k  numbers,  i.e.  7,8,  and  9  from  the  last  k  =  3  arrangements.  The 

deletion  of  certain  symbols  from  the  set  1,  2,  .  .  . ,  v  is  known  as  "variety  cutting". 
2 

For  k  =25  and  k  =  5  partially  balanced  incomplete  block  designs  may  be 

*  + 

constructed  for  v  =  10  and  k  =  2,  v  =  15  and  k  =3,  and  v  =  20  and  k'  =  4 
by  the  above  "variety  cutting"  procedure. 

2 

Also,  the  successive  diagonals  method  is  useful  for  v  =  k  in  incomplete 
blocks  of  size  k  for  any  odd  k  .  For  example,  for  v  =  225  and  k  =  15  four 
arrangements  or  squares  may  be  quickly  constructed  by  the  above  method.  Like¬ 
wise,  the  "variety  cutting"  procedure  may  be  utilized  to  obtain  2  or  3  arrange¬ 
ments  for  v  =  15p,  2  p  <_  15,  varieties  in  incomplete  blocks  of  size  p  . 

The  above  method  has  its  counterpart  in  constructing  mutually  orthogonal 
latin  squares  and  this  possibility  is  briefly  mentioned  in  Fisher  and  Yates  [1957] 
in  this  context.  Again  the  method  becomes  apparent  through  an  example.  First 
write  the  latin  square  in  standard  order  and  of  the  form  given  below  for  the  first 
square,  then  project  the  main  right  diagonal  of  the  preceding  square  into  the  first 
column  of  a  square,  and  then  write  the  symbol  order  in  the  same  manner  as  in  the 
first  square.  As  a  first  example,  let  the  order  of  the  latin  square  be  3;  the 
squares  are: 


first  square 


1 

2 

3 

2 

3 

1 

3 

1 

2 

second  square 


1 

2 

3 

3 

1 

2 

2 

3 

1 

third  square 


1 

2 

3 

1 

2 

3 

1 

2 

3 
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Thus,  the  main  right  diagonal  of  the  first  square  is  1,3,2  which  becomes  the 
first  column  of  the  second  square.  Then,  write  the  first  row  as  1,2,3,  the 
second  row  as  3,1,2,  and  the  third  row  as  2,3,1.  For  the  third  square, 
which  is  not  a  latin  square,  the  right  main  diagonal  of  the  second  square  is 
1,1,1  and  this  becomes  the  first  column  of  the  third  square}  the  rows  are  then 
completed.  If  we  then  take  the  right  main  diagonal  of  the  third  square,  we  ob¬ 
tain  the  first  square. 

As  a  second  illustrative  example,  the  five  squares  for  order  n  =  5  which 
are  constructed  by  successively  projecting  diagonals,  are: 

first  square  second  square  third  square 


The  fifth  square  is  not  a  latin  square  but  may  be  utilized  to  construct  the  first 
square  through  use  of  the  method  of  successive  projections  of  the  main  diagonals. 

The  method  may  be  utilized  for  any  odd  order  n  and  will  produce  q^-  1 
orthogonal  latin  squares  for  n  =  q^  q2  , . .  qg  where  q,  <  qi+1  and  qA  q2  ...  q 
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is  the  prime  power  decomposition  of  n  .  Thus,  for  n  =  15  =  3(5)  a  pair 
(q^-1  =  3-1  =  2)  of  orthogonal  latin  squares  is  easily  produced.  For  n  =  35  = 
5(7),  a  quartet  of  mutually  orthogonal  latin  squares  is  readily  produced  by  the 
projecting  diagonals  method. 


#1030 
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VII.  Relations  Between  Complete  Confounding  and  Simple  Qrthomorphlsms 


We  shall  illustrate  the  ideas  by  going  through  a  complete  example  taking 
2 

n  =  12  =  2  X  3  .  For  this  purpose  we  take  the  ring  of  12  elements  (obtained 
by  utilizing  Raktoe's  (1969]  results)  as  follows: 


GF(2Z)  _^3 

0  0 
1  =*  3 

x  3x 

x+1  3x+3 


GF(3)  ]± 

0  0 

1  =*  4 

2  2 


Rl2  =  I3  ®  I4  s  {0,1,  2,  3,4,5, 3x,  3x+l,  3x+2,  3x+3,  3x+4,  3x+5} 

R^  is  a  commutative  ring  under  addition  and  multiplication  (mod  (6, 

2 

3x  +  3x  +  3)  in  the  following  sense: 

e.g.  :  (a).  (3x+3)  +  (3x+4)  =  6x  +  ?  =  1;  here  we  have  to  reduce 

only  mod  6  to  get  the  answer. 

(b).  (3x+l)  •  ( 3x+  4 )  =  9x2  +  15x  +  4 

a  3x2  4  3x+  4  =  1  . 
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Explicitly,  to  facilitate  arithmetic,  the  addition  and  mutliplication  of  these  12 


elements  ares 


+ 

0 

1 

2 

3 

4 

5 

3x 

3x+l 

3x42 

3x+3 

3x44 

3x4  5 

0 

0 

1 

2 

3 

4 

5 

3x 

3x4 1 

3x+2 

3x4  3 

3x+4 

3x4  5 

1 

2 

3 

4 

5 

0 

3x4  1 

3x42 

3x+3 

3x+4 

3x+5 

3x 

2 

4 

5 

0 

1 

3x+2 

3X+3 

3x+4 

3x+5 

3x 

3x4 1 

3 

0 

1 

2 

3x+3 

3x+4 

3x45 

3x 

3x+l 

3x+2 

4 

2 

3 

3x+4 

3x+5 

3x 

3x+ 1 

3x+2 

3x4  3 

5 

4 

3x4  5 

3x 

3x+l 

3x42 

3x+3 

3X4  4 

3x 

0 

1 

2 

3 

4 

5 

3x+l 

2 

3 

4 

5 

0 

3x+2 

4 

5 

0 

1 

3x+3 

0 

1 

2 

3x+4 

2 

3 

3x+5 

4 
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Now  associate  with  a  latin  square  of  order  12  the 


32X42  =  [3X4]  X  [3X  4] 


=  12  X  12  lattice  square  with  the  following  breakdown  of  the  143  degrees  of 
freedom  corresponding  to  a  four-factor  factorial: 


A4 

2 

c3 

3 

B4 

2 

D3 

3 

a4b4 

2 

c3d3 

3 

aV 

2 

„3  3x 

C  D 

3 

_3_ 3x+3 

C  D 

3  j 

r 

,  4  3 

A  C 

6 

4  3 

B  C 

6 

4  4  3 

ABC 

6 

,4  2  3 

ABC 

6 

.43 

A  D 

6 

4  3 

B  D 

6 

,4  4  3 

A  B  D 

6 

aVd3 

6 

.433 

A  C  D 

6 

4  3  3 

B  CD 

6 

aVc3d3 

6 

a4b2c3d3 

6 

.  4  3_  3x 

A  C  D 

6 

-.4  3  3x 
BCD 

6 

aVc3d3x 

6 

a4b2c3d3x 

6 

a4c3d3x+3 

6 

_4  3_3x+3 
BCD 

6 

a4b4c3d3x+3 

6 

a4b2c3d3x+3 

6 

For  any  row  or  column  confounding  we  need  to  confound  effects  totaling  up  to  11 
degrees  of  freedom.  There  are  natural  candidates  available.  In  fact,  we  may 
choose  for  our  first  lattice  square  the  confounding  scheme  in  many  ways.  A 
scheme  resulting  in  a  pair  of  orthogonal  latin  squares  is  the  following: 
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Confounding  scheme  of  our  12  X  12  lattice  square 


#1030 
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LATIN  SQUARE  2;  Treatments  identified  with  ABC 


Using  the  complete  confounding  approach  as  outlined  above,  one  can  construct 
2 

min  [(2  -1),  (3-1)]  =  2  mutually  orthogonal  latin  squares  and  no  more  as  can 
easily  be  observed  from  the  degrees  of  freedom  table. 

From  the  multiplication  table  of  our  ring  R^,  we  observe  that  1,  5,  3x+l, 
3x+2,  3x+4  and  3x+5  are  the  6  non-zero  divisors  (i. e.  elements  with  multiplica¬ 
tive  inverses).  Following  Bose,  Chakravarti  and  Knuth  [I960]  ,  we  consider 

simple  automorphisms  in  of  the  form: 

* 

a(r)  =  r  r 

where  r  is  a  given  fixed  element  having  a  multiplicative  inverse  (because  only 
these  elements  are  capable  of  producing  automorphisms  of  R^)  •  Let  now  our 
aim  be  to  produce  two  orthomorphisms  which  in  turn  will  produce  an  0(12,  2) 
set.  For  this  purpose  consider  the  automorphisms: 

I(r)  =  r 

* 

a(r)  =  r  •  r  . 

Now  I  is  orthogonal  to  a  which  implies  the  condition  that  r  in  the  equation 

>!< 

[r  •  r-r]  =  c  has  a  unique  solution  for  every  c  of  R^  (see  Mann  [1949], 
pp.  103-105).  In  our  setting  this  means  that  r(r  -  1)  =  c  has  a  unique  solution, 

)j<  —  J 

i.  e. ,  r(r  +  5)  =  c  has  a  unique  solution  which  in  turn  implies  that  (r  +  5) 
exists  in  R  . 

1m 

j]c 

Substituting  in  the  values  of  r  we  see  that: 

[1  +  5]  *  does  not  exist  in  R^ 

[5  +  5]  ~ ^  "  "  "  "  " 

( ( 3x+l]  +  5]  *  does  not  exist  in  R^ 
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[  ( 3x4  2 )  -4  5]  *  =  [3x41]  '  'xists  m  R 

[(3x44)  4  5]  *  does  not  exist  in  R 

[(3x4  5  )  4  5]  -  [3x44j  ^  exists  in  R  . 

Hence  we  have  obtained  two  pairs  of  orthomorphisms  nar  -  1  : 

I(r)  =  r  l(r 

and 

a  (r )  =  ( 3x4  2  )r  a  (r )  -  t  4x )  5  )r  . 

J.  « 

The  0(12,  2)  set  presented  ibove  using  complete  confounding  corresponds  to 
the  first  pair  of  maps.  It  may  bo  asily  shown  that  sir. pie  m  g  >t  the  typ" 

j, 

Q'(r)  =  r  •  r  lead  to  0(12,  2)  sets  or  in  general  to  an  On.  a)  t.  when 
n.  n_  n,  k  n, 

a  =  min  (p.  -1,  p  -1,  . . .  ,  p  - 1)  and  n  =  |  |  p.  so  that  pj  tc  con- 

1  d  K  1=1  J 

founding  approach  is  equivalent  to  the  construction  of  a  set  of  a  orthomorphism:  . 
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VIII.  Some  Remarks  on  “Orthomorphism 11  Construction  of  Q(n,  t)  Sets 


In  1959  and/or  I960  E.  T.  Parker  showed  by  a  combination  of  classifica¬ 
tion  of  cases  (with  considerable  elimination  of  isomorphic  repetition  accessible 
to  cut  down  on  computer  time,  followed  by  computer  runs)  that  there  are  obtain¬ 
able  only  5  orthogonal  latin  squares  of  order  12,  all  restricted  to  be  copies  of 
the  non-cyclic  Abelian  group  with  latin  squares  related  by  row  permutations. 

(Some  researchers  [6,  31]  call  this  the  method  of  orthomorphisms.  Parker  considers 
this  no  method,  but  only  based  on  freaks  of  luck;  further,  Parker  feels  that 
"orthomorphism"  admits  no  precise  definition. ) 

Parker  made  another  finding,  also  by  hard  classification  of  cases  followed 
by  computer  runs,  which  Marshall  Hall  feels  is  more  important  than  that  cited 
above.  No  pair  of  order- 12  orthogonal  squares  of  the  type  mentioned  can  be  ex¬ 
tended  to  a  complete  set  of  any  sort;  i.e.  ,  further  orthogonal  latin  squares  are 
allowed  to  be  completely  general. 

What  might  be  obtainable  for  orthogonal  squares  of  order  20  in  like 
fashion,  row-permuting  the  non-cyclic  Abelian  group  of  order  20,  is  an  interest¬ 
ing  matter  for  speculation  —  conceivably  one  might  even  produce  a  complete  set 
(19  orthogonal  squares  equivalent  to  a  plane).  Knuth  and  Parker  discussed  the 
problem  about  1963,  and  concluded  that  exhaustive  search  is  out  of  the  question; 
still  a  fortunate  sample  of  cases  might  produce  an  attractive  result. 

In  1960  Parker  looked  at  the  row-permutation  ("orthomorphism"  of  Bose 
and  Mendelsohn)  approach  for  the  group  of  order  15,  and  proved  by  Hasse- 
Minkowski  invariants  that  a  complete  set  could  not  be  so  obtained.  He  dropped 
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further  work;  but  some  persistence  could  quite  possibly  yield  as  much  as  five 
orthogonal  squares  of  order  15  . 

A  hybrid  attack  on  order  15  or  20  might  be  undertaken  by  an  ambitious 
investigator.  (The  facts  for  order  12  mentioned  above  rule  out  chances  here.  ) 
One  might  produce  sets  of  orthogonal  latin  squares  of  row-permuted  group  type, 
using  automorphisms  of  the  group  latin  squares  to  eliminate  —  or,  that  failing, 
reduce  —  isomorphic  repetition.  It  would  not  be  shrewd  to  program  a  computer 
to  produce  all  transversals  of  a  group  latin  square,  for  running  time  and  output 
would  be  excessive;  then  for  any  hint  of  efficiency  it  would  be  necessary  to  turn 
about  and  do  a  reduction  on  the  computer  output.  After  a  set  of  row-permuted 
latin  squares  (possibly  exhaustive  for  order  15,  but  almost  certainly  only  a 
sample  for  order  20  )  large  enough  that  computer  searching  would  require 
realistic  amounts  of  time,  one  might  proceed  with  the  next  step.  Produce  all 
transversals  of  the  set  of  orthogonal  squares  by  computer,  then  fit  these  together 
in  all  possible  ways  (again  by  computer)  to  form  orthogonal  mates  of  the  preceding 
set  of  orthogonal  latin  squares.  Unlike  Parker's  assertion  above  about  complete 
sets  of  order- 15  squares,  there  is  no  known  argument  implying  impossibility  of 
producing  14  orthogonal  latin  squares  of  order  15  by  this  hybrid  attack. 
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IX.  Oval  Construction  of  Q(n,  t)  Sets 


The  approach  to  construction  of  finite  projective  planes  used  here  differs 

from  the  ones  known  in  the  literature.  The  mam  idea  is  to  make  use  of  the 

maximum  number  of  points  no  three  on  one  line  in  a  finite  projective  plane  of 

even  order  called  henceforth  an  oval.  The  oval  cannot  consist  of  more  than 

n  +  2  points  in  a  plane  of  order  n  .  This  is  obvious  since  through  each  point 
pass  n  +  1  lines  and  the  lines  through  any  point  of  the  oval  can  contain  at  most 
one  more  point  of  the  oval.  On  the  other  hand,  if  a  plane  of  order  n  does  include 
an  oval  consisting  of  the  maximum  number  of  points,  namely  n  +  2,  then  the 
lines  of  the  plane  can  be  classified  into  two  categories  in  respect  to  this  oval. 
One  category  consists  of  lines  intersecting  the  oval  in  two  points  called  secants, 

the  other  of  lines  having  no  points  of  the  oval  called  non- intersectors.  The 

number  of  secants  is  clearly  (n+2^(n+l)  ancj  number  of  non- inters ectors 

is  .  Through  each  of  the  n2-l  points  which  do  not  belong  to  the 

oval  there  would  have  to  pass  secants  and  ^  non- inters  ectors. 

Hence,  n  must  be  even. 
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It  is  well-known  that  removing  one  line  from  the  plane,  usually  called 

2 

the  line  at  infinity,  the  remaining  n  +  n  lines  can  be  arranged  into  2n  lines 

passing  through  two  points  at  infinity  which  are  arbitrary  up  to  notation  and 

2 

coordinatization  of  the  plane,  and  n  -n  lines  belonging  to  n-1  mutually 
orthogonal  latin  squares.  If  the  line  at  infinity  is  chosen  to  be  a  secant  and 
there  are  2n  lines,  the  lines  pass  through  the  two  points  of  the  oval  such 

that  each  of  the  n-1  latin  squares  consists  of  ^  secants  and  ^  non¬ 
intersectors  passing  through  each  of  the  n-1  points  at  infinity  other  than  the 

points  of  the  oval.  The  2n  lines  correspond  to  the  rows  and  columns  of  the  latin 
square. 

Using  the  described  method,  it  v/as  assumed  that  a  plane  of  order  10 
exists.  Under  this  assumption  21  lines  could  be  exhibited  arbitrarily  up  to 

notation.  Out  of  these  lines  one  was  taken  to  be  the  line  at  infinity  and  the 
remaining  20  used  to  coordinatize  the  plane.  Then  by  trial  and  error  twenty 
more  lines  were  found  which  formed  two  orthogonal  latin  squares.  The  method 
used  to  construct  these  squares  differs  from  the  one  described  in  the  literature. 
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Unfortunately  no  more  squares  could  be  found  using  this  method  and  a  high  speed 


computer  established  that  the  two  squares  did  not  yield  an  additional  mutually 
orthogonal  mate.  Clearly  it  could  happen  that  the  choice  of  the  first  two  was 
unfortunate.  The  same  method  was  applied  to  the  plane  of  order  12  .  Here 

the  trial  and  error  method  failed  to  produce  even  two  orthogonal  squares.  It 
may  be  worthwhile  to  remark  that  the  construction  of  the  plane  and  consequently 
the  search  for  orthogonal  latin  squares  does  not  require  the  assumption  that 
the  oval  consists  of  the  maximum  number  of  points  n  +  2  .  However,  if  the 
plane  does  not  include  an  oval  consisting  of  n  +  2  points  the  lines  could  not 
be  classified  into  two  categories  only  and  this  complicates  the  construction  of 

the  plane.  Let  us  illustrate  the  method  in  the  case  n  equals  10  .  It  is 
easy  to  show  that  in  this  case  the  oval  must  consist  of  at  least  6  points. 
However,  the  case  of  an  oval  of  6  points  would  be  ignored  since  in  this 
case  every  quadrangle  would  have  to  have  collinear  diagonals.  On  the  other 
hand,  a  plane  of  order  ten  must  be  non-Desarguessian  and  hence  must  contain 
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a  nondegenerate  quadrangle  with  noncollinear  diagonals.  Suppose  that  the 

plane  contains  a  quadrangle  with  noncollinear  diagonals  and  suppose  that  the 

plane  contains  an  oval  consisting  of  seven  points  then  the  104  points  of  the 

plane  which  do  not  belong  to  the  oval  could  be  classified  into  three  categories: 

(i)  points  lying  on  3  secants,  1  tangent  7  nonintersecfors 

(ii)  "  "  "  2  3  6 

(iii)  "  "  "  1  5  5 

Let  us  name  the  number  of  points  in  each  category  by  x,  y,  z  respectively. 
Clearly  x  +  y  +  z  =  104  . 

Counting  the  intersections  of  the  secants  and  the  tangents  we  get  the 
further  equations: 

3x  +  y  =  105 
3y  +  lOz  =  525 

The  unique  solutions  of  this  system  of  equations  are  x  =  20,  y  =  45,  z  =  39  . 
One  could  start  the  construction  of  the  plane  under  the  present  assumption  and 

investigate  the  possibilities  of  obtaining  orthogonal  latin  squares  in  this  way. 
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X.  Code  Construction  of  Q(n,  t)  Sets 


Given  an  n-symbol  alphabet,  e.g.  ,  l,2,...,n,  and  a  set  of  k-tuples 

of  the  n  symbols,  we  denote  the  set  of  all  k-tuples  by  C,  .  This  set  may 

k  ^  n 

be  thought  of  as  a  vector  space  or  as  a  k-dimensional  hypercube  with  edges  of 

length  n  .  Any  subset  of  C,  is  denoted  a®  a  block  code  with  a  block  length 

\  j  n 

of  k  .  The  elements  of  the  subset  are  denoted  as  code  words.  The  number  of 
symbols  by  which  any  two  code  words  differ  is  called  the  Hamming  distance. 

If  any  pair  of  code  words  in  the  subset  differs  by  a  Hamming  distance  of  at 
least  r,  the  block  code  is  called  a  distance  r  code .  A  distance  r  code  is 
called  an  (r-l)/2-error  correcting  code  because  fewer  than  (r-l)/2  changes 
leaves  the  word  closer  to  its  original  form  than  to  any  other  code  word  in  the 
subset.  For  similar  reasons,  a  distance  r  code  has  also  been  designated  as 
an  (r-l)-error-dectecting  code. 

In  an  interesting  paper,  Golomb  and  Posner  [1964]  discuss  the  relation- 

2 

ships  between  a  subset  of  n  code  words  and  an  0(n,t)  set  and  relate  these 

2 

to  ideas  developed  from  a  consideration  of  a  set  of  n  super  rooks  of  power  t 
t+2 

o-j  th  v  n  chessboard  such  that  no  two  super  rooks  attack  each  other.  The 

tus.i  concepts  of  rook  domains  and  rook  packing  were  found  to  be  very  useful  in 

providing  a  geometrical  view  of  the  results. 

2 

Any  subset  of  n  words  from  C,  which  forms  a  single-error-detecting 

3,n 

code  may  be  used  to  construct  a  latin  square  of  order  n  as  any  pair  of  the 

2 

triples  differs  by  at  least  two  symbols.  Likewise,  any  subset  of  n  words 
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from  C^+2  n  with  a  Hamming  distance  of  t  +  1  may  be  utilized  to  construct 
an  0<n,t)  set.  These  results  are  embodied  in  the  following  theorem  (from 
Golomb  and  Posner  [1964] ): 

Theorem  X.l.  The  following  three  concepts  are  equivalent: 

i)  an  0{n,t)  set. 

2  t+2 

ii)  A  set  of  n  nonattacking  super  rooks  of  power  t  on  the  n  board. 

2  / 

For  even  t,  also  the  following,  a  set  of  n  super  rooks  of  power  t/2  on  the 
t+2 

n  board  such  that  no  cell  is  attacked  twice;  that  is,  such  that  the  rook 
domains  are  nonoverlapping. 

2 

iii)  A  distance  t  +  1  code  of  block  length  t  +  2  with  n  words  from  an 
n-symbol  alphabet. 

For  those  interested  in  code  construction,  reference  may  be  made  to 

Mann  [1968]  and  Peterson  [1961]  and  the  literature  citations  therein.  We  shall 

2 

merely  illustrate  the  method  of  construction  of  an  0(n,t)  set  from  n  words 

of  length  t  +  2  and  Hamming  distance  t  +  1  through  an  example.  Let  n  =  3 

2 

and  t  =  2  .  Then  the  n  =9  code  words  with  length  4  and  Hamming  distance 
3  and  the  corresponding  la  tin  squares  are: 

latin  squares  of  order  3 


0  12  0  12 


0000 

0111 

0222 

0 

0 

1 

2 

0 

0 

1 

2 

1012 

1120 

1201 

to  produce 

1 

1 

2 

0 

and 

1 

2 

0 

1 

2021 

2102 

2210 

2 

2 

0 

1 

2 

1 

2 

0 

where  the  first  symbol  corresponds  to  row  number,  the  second  to  column  number, 
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the  third  to  symbols  in  the  first  la  tin  square,  and  the  fourth  to  symbols  in  the 
second  latin  square.  The  two  latin  squares  form  an  0(3,2)  set.  Note  that  any 
pair  of  the  quadruples  differs  in  at  least  three  symbols. 

The  analogy  of  the  above  with  many  of  the  concepts  from  fractional  repli¬ 
cation  and  orthogonal  arrays  is  immediately  apparent.  The  equivalences  of  many 
of  the  results  in  these  fields  need  to  be  systematically  noted  much  in  the  same 

manner  that  Golomb  and  Posner  [1964]  note  various  equivalences  among  0(n,t) 

2  t+2 

sets,  error-correcting  codes ,  and  n  nonattac king  rooks  on  an  n  chess¬ 
board. 
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XI.  Pairwise  Balanced  Design  Construction  of  Q(n,  t)  Sets 

Central  to  the  constructions  of  orthogonal  latin  squares  of  Bose  and 

Shrikhande  [1959]  and  of  Parker  [1959,  I960]  is  the  following  which  might  be 

called  a  "Folk  theorem,  "  being  credited  to  no  specific  investigator:  From  a  set 

2 

of  t  orthogonal  latin  squares  of  order  n  one  may  produce  a  set  of  n  ordered 
(t+2)-tuples  on  n  symbols  such  that  each  pair  of  distinct  positions  contains 
each  ordered  pair  of  symbols  (exactly  once);  the  converse  construction  can  also 
be  carried  out.  (Some,  such  as  Bose,  prefer  to  call  the  set  of  ( t+2)-tuples  an 
orthogonal  array.  )  There  is  nothing  difficult  to  prove  in  this  construction.  Two 
arbitrary  positions  in  the  (t+2)- tuples  are  identified  with  row  and  column  indices 
in  matrices,  and  each  other  position  with  entries  in  one  of  ti.o  matrices.  The 
equivalence  between  orthogonality  of  latin  squares  and  the  conditions  on  the 
(t+2)-tuples  is  then  fairly  apparent. 

Parker  [I960]  contributed  the  following  to  the  construction  of  orthogonal 


Let  the  3m  +  1  symbols  be  X^,. . .  ,Xm  and  the  residue  classes  (mod  2m  4  1). 
Form  the  (latin  square)  array 

Xi  0  i  -i 

0  X,  -i  i 

A:  1 

i  -i  Xi  0 

-i  i  0  XA 
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For  each  i,  1  :Li<_m,  each  row  of  A  an  ordered  quadruple.  In  turn,  the 
list  of  quadruples  is  built  up  by  adding  each  integer  (mod  2m  +  1)  to  all  four 
positions  at  once,  the  symbols  being  unchanged  by  the  addition.  The  set 
of  4m(2m  +  1)  ordered  quadruples  just  described  contains  in  each  pair  of  distinct 
positions  exactly  one  occurrence  of  each  ordered  pair  made  up  of  an  and  a 
residue  class  in  either  order,  and  of  each  ordered  pair  made  up  of  two  distinct 
residue  classes.  The  required  set  of  ordered  quadruples  is  completed  by  adjoin¬ 
ing:  i)  all  ordered  quadruples  (j,j,j,j),  j  =  0,  ...,  2m;  ii)  a  set  of  ordered 
quadruples  on  the  X,  symbols  corresponding  to  a  pair  of  orthogonal  latin  squares 
of  order  m  guaranteed  by  the  hypothesis  to  exist. 

Bose  and  Shrikhande  (1959,  published  1959  and  I960  partly  in  a  3-author 
paper  with  Parker)  developed  a  sequence  of  constructive  theorems  which  led  in 
steps  to  disproof  of  Euler's  conjecture  for  all  orders  4t  +  2  >  6  .  Their  central 
theorem  given  here  does  not  exhaust  their  methods,  but  virtually  all  their  results 
rest  on  this  theorem.  We  begin  with  a  definition.  A  pairwise  balanced  design. 

PB(n;  kj,. . .  ,k  )  is  a  collection  of  subsets  of  a  set  of  n  elements,  each  sub¬ 
set  having  number  of  elements  one  of  the  k^,  and  such  that  each  pair  of  distinct 
elements  in  the  set  of  n  occurs  in  a  unique  subset  of  the  PB  .  (Note:  unlike 
in  balanced  incomplete  block  designs,  the  subsets  of  a  PB  are  not  restricted  to 
have  equal  numbers  of  elements.  )  Now  for  the  main  theorem  of  Bose  and  Shrikhande. 
If  a  PB(n;  kj, . . .  ,k  )  exists,  and  for  each  i,  1  <_  i  <  t,  a  set  of  m  orthogonal 
latin  squares  of  order  k  exists .  then  a  set  of  m  -  1  orthogonal  latin  squares  of 
order  n  exists.  Loosely  speaking,  the  sets  of  ordered  tuples  for  each  subset 
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of  the  PB  are  constructed  and  these  fit  together  to  form  a  set  of  ordered  tuples 
for  the  full  set  of  n  elements.  The  decrease  from  m  to  m-1  orthogonal  latin 
squares  occurs  because  in  fitting  the  pieces  together  to  form  the  large  set  of 
ordered  tuples,  it  is  necessary  that  each  set  of  ordered  tuples  formed  from  a 
subset  of  the  PB  include  each  (i,i,...,i),  where  i  ranges  over  the  elements 
of  that  subset.  (It  is  sufficient  that  this  condition  be  fulfilled  in  the  construc¬ 
tion.  Thus  the  theorem  might  be  stated  in  slightly  stronger  form:  " If  . . .  1<  i<_t  , 
a  set  of  m  orthogonal  latin  squares  of  order  with  a  transversal  exists,  then 
a  set  of  m  orthogonal  latin  squares  of  order  n  exists.  ")  Now  for  a  more  nearly 

formal  version  of  the  proof.  If  there  exists  a  set  of  m  orthogonal  latin  squares 

2 

of  order  n,  then  there  exists  a  set  of  the  appropriate  sort  of  n  ordered 
(m-fl)-tuples  with  each  symbol  repeated  in  an  (m+l)-tuple  m  +  1  times.  (The 
condition  mentioned  is  satisfied  with  (m±2)- tuples  if  the  set  of  orthogonal  latin 
squares  has  a  transversal. )  One  need  simply  put  together  the  ordered  tuples  on 
each  subset  of  the  PB  in  turn,  subject  to  the  important  condition  that  within 
each  subset  of  the  PB,  each  tuple  of  repetitions  of  each  symbol  be  included. 
Carrying  this  out  on  the  alphabet  of  the  symbols  in  each  subset  of  the  PB,  one 
has  the  construction  for  the  set  of  orthogonal  latin  squares  in  the  conclusion; 
each  ordered  tuple  of  a  repeated  symbol  among  the  n  is  used  only  once. 

A  representative  and  very  interesting  example  (Bose  and  Shrikhande  in¬ 
formed  Parker  that  this  was  the  first  case  of  disproof  of  Euler's  conjecture  pro¬ 
duced  in  their  joint  work  at  a  blackboard )  yields  5  mutually  orthogonal  latin 
squares  of  order  50  via  the  PB  construction.  One  forms  the  affine  plane  of 
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order  7,  then  adjoins  exactly  one  ideal  point  on  each  line  of  one  class  of  paral¬ 
lel  lines.  This  yields  a  PB(50 ;  8,7).  Since  there  exist  6  orthogonal  latin 
squares  of  each  order  8  and  7,  there  exist  6-1=5  orthogonal  latin  squares 
of  order  50  . 

There  is  a  limitation  on  the  Bose-Shrikhande  PB  construction.  Aside  from 
trivial  PB  designs,  having  a  single  subset  of  all  elements,  any  PB  has  a  sub¬ 
set  with  at  most  one  more  element  than  the  square  root  of  the  number  of  elements 
in  the  large  set.  Thus  other  techniques  are  requisite  to  produce  more  than  si n 
orthogonal  latin  squares  of  order  not  a  prime-power. 


-78- 


#1030 


XII.  Product  Composition  of  Q(n,  t)  Sets 


About  70  years  ago,  for  the  first  time,  Tarry  [1899}  in  his  half-page  note 
asserted  that  if  there  exists  an  0{a,  2)  set  and  if  there  exists  an  0(b,2)  set 
then  there  exists  an  0(ab,2)  set.  He  exhibited  the  following  0(12,2)  set,  by 
composing  two  0(3,2)  and  0(4,2)  sets,  to  demonstrate  the  truth  of  his  asser¬ 
tion.  Note  that  in  the  following  square  the  set  of  first  integers  belong  to  one 
latin  square  and  the  set  of  second  integers  belong  to  the  second  latin  square. 

No  more  description  is  given  by  Tarry. 
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Tarry  did  not  observe  any  generalization  of  his  method.  Perhaps  this  was 
due  to  the  fact  that  he,  like  so  many  other  researchers,  was  only  concerned  with 
sets  of  type  0(n,  2)  .  Probably  he  was  not  aware  of  the  existence  of  a  larger  set. 
About  23  years  later  MacNeish  [1922]  demonstrated: 
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1)  The  existence  and  a  construction  of  an  0(n,  n-1)  set  for  n  a  prime  or 
prime  power  integer. 

2)  A  generalization  of  Tarry's  procedure  viz. ,  if  there  exists  an  0(a,r)  set  and 
if  there  exists  an  0(b,r)  set  then  there  exists  an  0(ab,  r)  set. 

a.  a  a 

12  t 

3)  By  a  successive  application  of  1)  and  2)  he  showed  that  if  n  =  p^  p2  . . .  Pt 
is  the  prime-power  decomposition  of  n  then  there  exists  an  0(n,  r)  set  where 

“i 

r  =  min{p.  -  1,  i  =  1, 2, . . .  ,t  }  . 

MacNeish  could  not  embed  his  0(n,  r)  set  generated  in  3)  in  a  larger 
set.  This  unsuccessful  attempt,  reinforced  by  Filer's  conjecture,  led  MacNeish 
to  prove  (erroneously)  geometrically  that  0(n,z)  sets  do  not  exist  for  z  >  r  , 
and  therefore,  as  a  confirmation  of  Euler's  conjecture.  The  preceding  argument 
of  MacNeish  is  known  as  MacNeish' s  conjecture  in  the  literature.  By  constructing 
an  0(21,3)  set  Parker  [1959]  gave  a  counter  example  to  MacNeish's  conjecture. 
Later  Bose,  Shrikhande,  and  Parker  [1960]  completely  demolished  Euler's  con¬ 
jecture  except  for  n  =  6  .  It  should  be  mentioned  that  MacNeish's  conjecture 
has  not  been  totally  disproved  yet.  For  instance,  no  one  as  yet  as  far  as  we  know, 
has  constructed  an  0(15,4)  set  (an  0(15,3)  set  is  given  in  section  V  for  the 
first  time)  or  an  0(20,3)  set.  We  believe  that  MacNeish  should  be  given  sub¬ 
stantial  credit  for  his  non-erroneous  contributions.  It  is  to  be  regretted  that 
MacNeish  is  often  cited  in  the  literature  only  for  his  false  conjecture. 

Even  though  Tarry  and  MacNeish  did  not  attach  any  name  to  their  pro¬ 
cedure,  it  is  not  difficult  to  see  that  it  is  the  method  of  Kronecker  product  of 
matrices.  Therefore,  we  can  state,  more  formally,  their  results  as  follows: 
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Theorem  (Tarry- MacNelsh).  If.  {A^,  A., . A  }  is  an  0(n,r)  set  and  if 

(Br  B2 . Br)  is  an  O(m.r)  set,  then  {Aj® Bj,  Ag  ®  B2>  . . .  ,  Ar  ®  Br }  , 

where®  denotes  the  Kronecker  product  operation  of  matrices,  is  an  0(nm,r) 
set. 

The  preceding  arguments  clearly  support  the  choice  of  the  title  for  this 
section  and  is  in  contrast  to  the  choice  of  the  name  for  the  procedure  given  in 
section  XIII. 
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XIII.  Sum  Composition  Construction  of  Q(n,  t)  Sets 
XIII.  1.  Introduction 

Perhaps  one  of  the  most  useful  techniques  for  the  con¬ 
struction  of  combinatorial  systems  is  the  method  of  composition.  To  mention 
some,  here  are  few  well-known  examples:  1)  If  there  exists  a  set  of  t  ortho¬ 
gonal  latin  squares  of  order  n^  and  if  there  exists  a  set  of  t  orthogonal  latin 
squares  of  order  n  ,  then  there  exists  a  set  of  t  orthogonal  latin 
squares  of  order  n^  .  2)  If  there  are  Steiner  triple  systems  of  order  and 
v2,  there  is  a  Steiner  triple  system  of  order  v  =  .  3)  If  and  H2  are 

two  Hadamard  matrices  of  order  n^  and  n2  respectively,  then  the  Kronecker 
product  of  Hj  and  is  a  Hadamard  matrix  of  order  n^2  .  4)  If  Room 
squares  of  order  n^  and  n2  exist,  then  a  Room  square  of  order  n^2  exists. 

5)  If  BIB  (v  ,  k,  \j)  and  BIB  (v2, k,  exist  and  if  k,  then 

2 

BIB  (vjv2>  k,\  2>  exists  where  f(\  2v2)  denotes  the  maximum  number  of  con- 

2 

straints  which  are  possible  in  an  orthogonal  array  of  size  \  v  ,  with  v  levels, 

Ct  Ct  b 

strength  2,  and  index  \  .  6)  As  a  final  example,  the  existence  of  orthogonal 

w 

arrays  (  \  ^v*,  q^,  v.,  t),  i  =  1,  2,  . . . ,  r  implies  the  existence  of  the  orthogonal 
array  (Kv*,  q,v,t),  where  \  =  ^  Vs  vlV2**,Vr’  and  q=  min(qi>q2» 

•  • • >  dr)  • 

The  reader  will  note  that  each  of  the  above  examples  involved  a  product 
tvpe  composition.  The  method  that  we  will  describe  utilizes  a  sum  type  compo¬ 
sition,  by  means  of  which  one  can  possibly  construct  sets  of  orthogonal  latin 
squares  for  all  n  >  10  . 

XIII. 2.  Definitions 

In  the  sequel  by  an  0(n,  t)  set  we  mean  a  set  of  t  mutually  orthogonal 
latin  squares  of  order  n  . 
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a)  A  transversal  (directrix)  of  a  latin  square  L  of  order  n  on  an  n-set  2 


is  a  collection  of  n  cells  such  that  the  entries  of  these  cells  exhaust  the  set 
2  and  every  row  and  column  of  L  is  represented  in  this  collection.  Two  trans¬ 
versals  are  said  to  be  parallel  if  they  have  no  cell  in  common, 
b)  A  collection  of  n  cells  is  said  to  form  a  common  transversal  for  an  0(n,  t) 
set  if  the  collection  is  a  transversal  for  each  of  these  t  latin  squares. 

Two  common  transversals  are  said  to  be  parallel  if  they  have  no  cell  in  common. 
Example.  The  underlined  and  parenthesized  cells  form  two  parallel  common 
transversals  for  the  following  0(4,2)  set. 

r  1  2  (3)  4 

1(2)  l  4  3 

]  3  (4)  1  2 

v.  4  3  2  (1) 


1  2  (3)  4  ^ 

(4)  3  2  1 

2  (1)  4  3 

3  4  I  (2)> 


XIII.  3.  Composing  Two  Latin  Squares  of  Order  n^  and  n 2 

Avery  natural  question  in  the  theory  of  latin  squares  is  the  following: 
Given  two  latin  squares  and  L 2  of  order  n^  and  n2  (n^  >  n2)  respectively. 
In  how  many  ways  can  one  compose  and  L2  in  order  to  obtain  a  latin  square 
L3  of  order  m,  where  m  is  a  function  of  n^  and  n2  only?  This  question 
can  be  partially  answered  as  follows.  First,  it  is  well-known  that  the  Kronecker 
product  L  =  L.  ®  L  is  a  latin  square  of  order  m  =  n.n  irrespective  of  the 

jib  1  ^ 

combinatorial  structure  of  and  L2  .  Secondly,  we  show  that  if  has  a 
certain  combinatorial  structure,  then  one  can  construct  a  latin  square  L  of 
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order  n  =  n^  +  .  Naturally  enough  we  call  this  procedure  a  "method  of  sum 

composition". 

Even  though  our  method  of  sum  composition  does  not  work  for  all  pairs 
of  latin  squares,  it  has  an  immediate  application  in  the  construction  of  ortho¬ 
gonal  latin  squares  including  those  of  order  4t  +  2,  t  >  2  .  We  emphasize  that 
the  combinatorial  structure  of  orthogonal  latin  squares  constructed  by  the  method 
of  sum  composition  is  completely  different  from  those  of  known  orthogonal  latin 
squares  in  the  literature.  Therefore,  it  is  worthwhile  to  study  these  squares 
for  the  purpose  of  constructing  new  finite  projective  planes. 

We  shall  now  describe  the  method  of  "sum  composition".  Let  and 

be  two  latin  squares  of  order  n^  and  n2,  n^  >  n2  ,  on  two  non-intersecting 

sets  Sj  =  {aj,  a2,  . . .  ,  an  }  and  "L ^  =  {bpb2,  . . .  ,  b^  }  respectively.  If  Lj 

1  2 

has  n2  parallel  transversals  then  we  can  compose  with  L2  to  obtain  a 
latin  square  L  of  order  n  =  n^  +  n2  .  Note  that  for  any  pair  (n^,  n2),  there 
exists  Lj  and  L2  with  the  above  requirement,  except  for  (2,1),  (2,2),  (6,5) 
and  (6,6)  . 

To  produce  L  put  and  L2  in  the  upper  left  and  lower  right  corner 
respectively.  Call  the  resulting  square  C^,  which  looks  as  follows: 


Name  the  n2  transversals  of  in  any  manner  from  1  to  n2  .  Now  fill  the 
cell  (i,  n^  +  k),  k  =  l,2,...,n2,  with  that  element  of  transversal  k  which 
appears  in  row  i,  1*1,2,,..,^.  Fill  also  the  cell  (nj+k,)),  k  =  1, 2, . . .  ,n2  , 
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with  that  element  of  transversal  k  which  appears  in  column  j,  j  =  1, 2, . . .  ,n^  . 
Call  the  resulting  square  C,  .  Now  every  entry  of  C,  is  occupied  with  an 
element  either  from  2j  or  22,  but  C2  is  obviously  not  a  latin  square  on 
2j  U  22  .  However,  if  we  replace  each  of  the  n^  entries  of  transversal  k 
with  b^,  it  is  easily  verified  that  the  resulting  square  which  we  call  L  is  a 
latin  square  of  order  n  on  2^  U  T, ^  . 

The  procedure  described  for  filling  the  first  n^  entries  of  the  row  (column) 
nj  +  k  with  the  corresponding  entries  of  transversal  k  is,  naturally  enough, 
called  the  projection  of  transversal  k  on  the  first  n^  entries  of  row  (column) 

^  +  k  . 


We  shall  now  elucidate  the  above  procedure  via  an  example. 


Let  2  ^  = 


{1,2, 3,4,5},  2 2  =  {6,7,8}  , 


6  7  8 
and  L  =  786 

w 

8  6  7 


Note  that  the  cells  on  the  same  curve  in  form  a  transversal. 


1 

2 

3 

4 

5 

1 

2 

3 

4 

5 

1 

2 

3 

5 

1 

2 

3 

4 

5 

1 

2 

3 

4 

4 

5 

1 

4 

5 

1 

2 

3 

4 

5 

1 

2 

3 

2 

3 

4 

3 

4 

5 

1 

2 

and  C2 

3 

4 

5 

1 

2 

5 

1 

2 

2 

3 

4 

5 

1 

2 

3 

4 

5 

1 

3 

4 

5 

6 

7 

8 

1 

3 

5 

2 

4 

6 

7 

8 

7 

8 

6 

5 

2 

4 

1 

3 

7 

8 

6 

8 

6 

7 

4 

1 

3 

5 

2 

8 

6 

7 
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And  finally 


6 

7 

8 

4 

5 

1 

2 

3 

7 

8 

2 

3 

6 

4 

5 

1 

8 

5 

1 

6 

7 

2 

3 

4 

3 

4 

6 

7 

8 

5 

1 

2 

2 

6 

7 

8 

1 

3 

4 

5 

1 

3 

5 

2 

4 

6 

7 

8 

5 

2 

4 

1 

3 

7 

8 

6 

4 

1 

3 

5 

2 

8 

6 

7 

which  is  a  la  tin  square  of  order  8  on  2^  U  22  =  {1,2,. ..,8}. 

Remark.  Note  that  it  is  by  no  means  required  that  the  projection  of  transversals 

on  the  rows  and  columns  should  have  the  same  ordering.  Indeed,  for  the  fixed 

set  of  ordered  n^  transversals,  we  have  n2*.  choices  of  projections  on  columns 

and  n2'.  choices  of  projections  on  the  rows.  Hence  we  can  generate  at  least 
2 

(n2i  )  different  la  tin  squares  of  order  n  =  ni  +  n2  composing  and  L2  . 

XIII.  4.  Construction  of  0(n,2)  Sets  by  Method  of  Sum  Composition.  In  order 

In  order  to  construct  an  0(n,  2)  set  for  n  *  n^  +n2,  we  require  that  n^>  2n2  and 
there  should  exist  an  0(n2,2)  set,  and  an  0(n^,2)  set  with  2n2  parallel 
transversals.  It  is  easy  to  show  that  any  n  >  10  can  be  decomposed  in  at  least 
one  way  into  n^  +  n2  which  fulfill  the  above  requirements.  We  now  present  two 
theorems  which  state  that  for  certain  n  one  can  construct  an  0(n,2)  set  by  the 
method  of  sum  composition. 

Theorem  XIII.  4. 1.  Let  n^  =  p  >7  for  any  odd  prime  p  and  positive  integer 
a,  excluding  n^  =  13.  Then  there  exists  an  0(n,2)  set  which  can  be  constructed 
by  composition  of  0(n^,2)  and  0(n2,2)  sets  for  n2  =  ( n ^  —  1 )/ 2  and 
n  =  n^  +  n2  . 
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We  shall  first  give  the  method  of  construction  and  then  a  proof  that  the 
constructed  set  is  an  0(n,2)  set. 

Construction.  Let  B(r)  be  the  n^  X  n^  square  with  element  ra^+a.  in  its  (i,j) 
cell,  a afj ,  0  *  r  in  GF(n^),  i,  j  =  1,2,. . .  ,n^  .  Then  it  is  easy  to  see  that 
{ B(  1 ) ,  B(x),  B(y)},  y  =x  \  x  *  y,  is  an  0(^,3)  set.  Consider  the  n1  cells  in  B(l) 
with  +  a  =  k  a  fixed  element  in  GF(n^)  .  Then  the  corresponding  cells  in 
B(x)  and  B(y)  form  a  common  transversal  for  the  set  {B(x),  B(y)} .  Name  this 
common  transversal  by  k  .  It  is  then  obvious  that  two  common  transversals  k^ 
and  k^,  k^  #  k^  are  parallel  and  hence  {B(x),  B(y)}  has  n^  common  parallel 
transversals.  Now  let  {A^A^}  beany  0{n2,2)  set,  which  is  known  to  exist,  on 
a  set  fi  non-intersecting  with  GF(n^)  .  For  any  X  in  GF(n^)  we  can  find 
(n1*l)/2  pairs  of  distinct  elements  belonging  to  GFfn^)  such  that  the  sum  of 
the  two  elements  of  each  pair  is  equal  to  X.  Let  {S}  and  {T}  denote  the 
collection  of  the  first  and  the  second  elements  of  these  (n^-  1)/ 2  pairs  respec¬ 
tively.  Note  that  for  a  fixed  X  the  set  {S}  can  be  constructed  in  (n^-lHn^-  3) 

2 

...1  distinct  ways.  Now  fix  X  and  let  L  denote  any  of  the  (n  I  )  latin 
squares  that  can  be  generated  by  the  sum  composition  of  L(x)  and  A^  using 
transversals  determined  by  the  n^  elements  of  {S}  .  Let  L^  be  the  latin  square 
derived  from  the  composition  of  L(y)  and  A^  using  the  n^  transversals  de¬ 
termined  by  the  elements  of  {T}  and  the  following  projection  rule:  Project  trans¬ 
versals  tj,  i  -  1,2,...  ,n^  on  the  row  (column)  which  upon  superposition  of 
L^  on  L^  this  row  (column)  should  coincide  with  the  row  (column)  derived  from 
the  projection  of  the  transversal  X-t^  .  shortly  we  shall  prove  that  { 


forms  on  0(n,2)  set. 
#1030 
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The  preceding  arguments  shows  that  { }  can  be  constructed  non- 

2 

isomorphically  in  at  least  (n^-3)(n2!  )  [  n^(n^-l)(n^-3). . .  lj  ways.  For  instance 
in  the  case  of  n^  =  7,  there  is  at  least  12096  non-isomorphic  pairs  of  ortho¬ 
gonal  latin  squares  of  order  10  .  Therefore,  Euler  has  been  wrong  in  his  con¬ 
jecture  by  a  very  wide  margin. 

Note  that  we  can  construct  infinitely  many  pairs  of  orthogonal  latin  squares 
of  order  4t  +  2  by  the  method  of  theorem  XIII.  4. 1.  For  p  =  7  mod  8  and  a  odd 
p^  =  (8t  +  5)/ 3  .  Hence  n^  +  n2  =  4t  +  2  . 

Proof;  The  constructional  procedure  clearly  reveals  that: 

A.  1^  and  L2  are  latin  squares  of  order  n  on  GF(n^)  U  S3  . 

B.  Upon  superposition  of  on  the  following  are  true: 

b^.  Every  element  of  S3  appears  with  every  other  element  of  S3  . 
b2<  Every  element  of  S3  appears  with  every  element  of  GF(n^)  . 
b^.  Every  element  of  GF(n^)  appears  with  every  element  of  S3  . 

Therefore,  all  we  have  to  prove  is  that  every  element  of  GF(n^)  appears  with 
every  other  element  of  GF(n^)  .  To  prove  this  recall  that  B(x)  is  orthogonal  to 
B(y)  .  However,  since  we  removed  the  n2  transversals  from  B(x)  determined 
by  the  n2  elements  of  {S}  and  n2  transversals  from  B(y)  determined  by  the 
n2  elements  of  {T}  therefore  the  following  2n2n^  pairs  have  been  lost. 

(xo'i  +  a. ,  ya.  +  a, )  with  +  =  V  for  any  \  «  GF(n^),  \  ±  \  . 

We  claim  that  the  given  projection  rules  guarantee  the  capture  of  these  lost  pairs 
by  the  2n2n^  bordered  cells.  To  show  this  note  that  the  superposition  of  the 
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projected  transversal  s  from  B(x)  on  the  projected  transversal  t  =  X  -  s  from 
B(y)  will  capture  the  n^  pairs 

(xa,  +  Oj,  y a.  +  a^)  with  a.  +  a.  =  k  =  [y(X -s)  +  sJ/(l-»  y) 
if  these  transversals  have  beer  projected  on  row  border  and  the  n^  pairs 
(xa^  +  a, ,  ya^  +  a, )  with  a.  +  a,  =  k  =  [  s(y-l)  +  (s-  X  )(x—  1 ) J  /(y-x) 

if  these  transversals  have  been  projected  on  column  border.  Now  because 
k  +  k‘  =  X  and  if  s^  4  s^  then  k^  4  k^  and  kj  4  k^  hence  the  2n^n^  pairs 
which  have  been  resulted  from  the  projection  of  transversals  determined  by  {S} 
and  {T}  will  jointly  capture  the  2n^n^  lost  pairs  and  thus  a  proof. 

We  shall  now  clarify  the  above  constructional  procedure  by  an  example. 
Example.  Let  n^  =  7,  GF(7)  =  {0,1,2,.. .  ,6}  .  Then  for  x  =  2,  y  =  x  *  =  4 
we  have 


{ B(l) ,  B(2 ) ,  B(4) }  = 

0123456  0123456  0123456 

1234560  2345601  4560123 

2345601  4  5  60123  1234560 

3456012  6012345  5601234 

4560123  1234560  2345601 

5601234  3456012  6012345 

6012345  5601234  3456012 


For  n2  =  (n^lj/2  =3  let  S22  =  {7,8,9}  and 
7  8  9  7  8  9 

{A  , A  }=  8  9  7,  9  7  8  .  Finally  for  X  =  0,  {S}  =  { 1 , 2, 3}  and 

9  7  8  8  9  7 

{T}  =  {6,5,4}  we  have  {L^L^ 
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the  reader  can  easily  verify  that  {L  ,  L_}  is  an  0(10,2)  set. 

Remarks. 

1)  The  method  of  theorem  XIII.  4. 1  fails  for  n^  =  13  only  because  there  is  no 
0(6,2)  set.  Otherwise,  there  will  be  no  orthogonality  contradiction  on  the  other 
parts  of  Lj  and  with  their  6x6  lower  right  square  missing. 

2)  In  the  case  of  n^  =  7,  if  we  let  {S}  =  {0,1,3}  and  {T}  =  {2,4,5}  then  the 
requirement  y  =  x  *  is  not  necessary.  However  then  we  do  not  have  a  unified 
projection  rule  for  the  formation  of  L2  as  was  provided  for  the  case  y  =  x  *  by 
theorem  XIII.  4. 1.  To  give  the  complete  list  of  solutions  let  (a^,a2,a3)  and 
(b^b^jb^)  be  any  two  permutations  of  the  set  {8,9,10} .  If  we  project  trans¬ 
versals  (0,1.3)  on  the  rows  (a^a^a^)  and  columns  (b^,b2,b3)  in  the  forma¬ 
tion  of  Lj,  then  the  following  table  indicates  what  permutation  of  transversals 
{2,4,5}  should  be  projected  on  the  rows  (a^,a2,a3)  and  columns  (b^,b2,b3) 
in  the  formation  of  .  Obviously  these  permutations  will  be  a  function  of  the 
pair  (x,y)  . 


-90- 


#1030 


Pair 

(x,y) 

Rows 

Wa3 

Columns 

bl»b2,b3 

(2,3) 

4,  2,  5 

4,  2,  5 

(2,3) 

2,  5,  4 

2,  5,  4 

(2,4) 

2,  5,  4 

4,  2,  5 

(2,5) 

4,  2,  5 

4,  2,  5 

(2,6) 

2,  5,  4 

2,  5,  4 

(3,4) 

2,  5,  4 

2,  5,  4 

(3,5) 

2,  5,  4 

4,  2,  5 

(3,5) 

4,  2,  5 

5,  4,  2 

(3,5) 

4,  2,  5 

2,  5,  4 

(3,5) 

5,  4,  2 

2,  5,  4 

(3,6) 

4,  2,  5 

2,  5,  4 

(3,6) 

5,  4,  2 

4,  2,  5 

(4,5) 

2,  5,  4 

2,  5,  4 

(4,6) 

5,  4,  2 

4,  2,  5 

(4,6) 

2,  5,  4 

2,  5,  4 

(4,6) 

5,  4,  2 

5,  4,  2 

(This  table  is  by  no  means  exhaustive. ) 

The  reader  may  note  that  whenever  y  =  x  *  in  the  above  table  the  given 
solutions )  are  different  from  the  one  provided  by  the  method  of  theorem  XIII.  4. 1. 

Thus  we  can  conclude  that  any  pair  of  orthogonal  latin  squares  of  order  7 
based  on  the  GF(7,  can  be  composed  with  a  pair  of  orthogonal  latin  squares  of 


#1030 


-91- 


order  3  and  make  a  pair  of  orthogonal  latin  squares  of  order  10  .  In  addition, 
since  we  have  six  choices  for  (a^a^a^)  and  (b^,b2,b3)  hence  from  every 
line  in  the  above  table  we  can  produce  36  non-isomorphic  0(10,2)  sets  or 
16  X  36  =  576  sets  for  the  entire  table.  Since  all  these  pairs  are  non-isomorphic 
with  all  previous  pairs,  produced  by  theorem  XIII.  4.1,  thus  by  the  method  of  sum 
composition  one  can  at  least  produce  12,672  non-isomorphic  0(10,2)  sets. 

We  believe  that  for  other  values  of  n^  there  are  sets  of  {S}  and  {T} 
together  with  proper  projections  which  makes  the  restriction  y  =  x  *  unnecessary. 
Theorem  XIII.  4.  2.  Let  n^  =  2a  >  8  for  any  positive  integer  a  .  Then  there 
exists  an  0(n,2)  set  which  can  be  constructed  by  composition  of  0(n^,2) 
and  0(n2,2)  sets  for  n^  =  n^/2  and  n  =  n^  +  n^  . 

We  shall  here  give  only  the  method  of  construction.  A  similar  argument 
as  in  theorem  XIII.  4.1  will  show  that  the  constructed  set  is  an  0(n,2)  set. 
Construction.  In  a  similar  fashion  as  in  theorem  XIII.  4. 1  construct  the  set 
(B(l),  B(x),  B(y )}  over  GF(2a’)  .  Let  also  (ApA2}  beany  0(n2,2)  set,  which 

(X 

always  exists,  on  a  set  f2  non-intersecting  with  GF(2  )  .  For  any  X.  *  0  in 
GF(2W)  we  can  find  n^/2  pairs  of  distinct  elements  belonging  to  GF(2a)  such 
that  the  sum  of  the  two  elements  of  each  pair  is  equal  to  \  .  Let  {S}  and  {T} 
denote  the  collection  of  the  first  and  the  second  elements  of  these  n^/2  pairs 
respectively.  Note  that  for  a  fixed  \  the  set  {S}  can  be  constructed  in 
n^(n^-2)(n^-4). . .  1  distinct  ways.  Now  form  L^  from  the  sum  composition  of  B(x) 
and  Aj  and  L^  from  the  sum  composition  of  B(y)  and  using  the  same  pro¬ 
jection  rule  as  given  in  theorem  XIII.  4.1.  Now  {L^,L2}  is  an  0(n,2)  set. 
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Example.  Let  n  =  8,  GF(8)  =  {0,1,2,. ..  ,7}  with  the  following  addition  (+) 
and  multiplication  (X  )  tables: 


+ 

0 

1 

2 

3 

4 

5 

6 

7 

X 

0 

1 

2 

3 

4 

5 

6 

7 

0 

0 

1 

2 

3 

4 

5 

6 

7 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

6 

4 

3 

7 

2 

5 

1 

0 

1 

2 

3 

4 

5 

6 

7 

2 

2 

6 

0 

7 

5 

4 

1 

3 

2 

0 

2 

3 

4 

5 

6 

7 

J 

3 

3 

4 

7 

0 

1 

6 

5 

2 

3 

0 

3 

4 

5 

6 

7 

1 

2 

4 

4 

3 

5 

1 

0 

2 

7 

6 

4 

0 

4 

5 

6 

7 

1 

2 

3 

5 

5 

7 

4 

6 

2 

0 

3 

1 

5 

0 

5 

6 

7 

1 

2 

3 

4 

6 

6 

2 

1 

5 

7 

3 

0 

4 

6 

0 

6 

7 

1 

2 

3 

4 

5 

7 

7 

5 

3 

2 

6 

1 

4 

0 

7 

0 

7 

1 

2 

3 

4 

5 

6 

-1 

Then  for 

X 

2, 

y  : 

=  > 

= 

7 

we  have 

1B(1),  B(2), 

,  B(7 )}  - 

0 

1 

2 

3 

4 

5 

6 

7 

0  1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

2 

;  4 

5 

6 

7 

1 

0 

6 

4 

3 

7 

2 

5 

2  6 

0 

7 

5 

4 

1 

3 

7 

5 

3 

2 

:  6 

1 

4 

0 

2 

6 

0 

7 

5 

4 

1 

3 

3  4 

7 

0 

1 

6 

5 

2 

1 

0 

6 

4  3 

7 

2 

5 

3 

4 

7 

0 

1 

6 

5 

2 

4  3 

5 

1 

0 

2 

7 

6 

2 

6 

0 

1 

'  5 

4 

1 

3 

4 

3 

5 

1 

0 

2 

7 

6 

5  7 

4 

6 

2 

0 

3 

1 

3 

4 

7 

0  1 

6 

5 

2 

5 

7 

4 

6 

2 

0 

3 

1 

6  2 

1 

5 

7 

3 

0 

4 

4 

3 

5 

] 

0 

2 

7 

6 

6 

2 

1 

5 

7 

3 

0 

4 

7  5 

3 

2 

6 

1 

4 

0 

5 

7 

4 

6  2 

0 

3 

1 

7 

5 

3 

2 

6 

1 

4 

0 

1  0 

6 

4 

3 

7 

2 

5 

6 

2 

1 

C 

i  7 

3 

0 

4 

For  n^  =  n^/2 


4  iet  S2  =  {A,  B,  C,  D}  and 


A  B  C  D 

{A  A  }  =  B  A  D  C 
1  F  2'  C  D  A  B 
D  C  B  A 


A  B  C  D 
D  C  B  A 
B  A  D  C 
C  D  A  B 


Finally  for  \  =  5,  {S}  =  {0,1, 3, 4}  and  {T}  =  {5, 7,6,2}  we  have  {L^H 


A  B  2  C  D  5  6  7 
BA0DC413 
34A01DBC 
C  D  5  A  B  2  7  6 
DC4BA031 
62D57ACB 
75B26CAD 
10C43BDA 

0134  01D34ACB 

6257  75C26BDA 

7526  DC6BA725 

1043  26B75CAD 

2675  34A01DBC 

3401  AB5CD276 

4310  CD4AB031 

5762  B  A  1  D  C  3  0  4 

5  7  6  2 

0  13  4 

3  4  0  1 
10  4  3 

7  5  2  6 

4  3  10 

6  2  5  7 

2  6  7  5 

0  6  7  1  2  3  4  5 

ABCD  4  2  >  6  3  5  1  0 

ABCD 

2  1  3  6  0  7  5  4 

BADC  63042157 

DCBA 

4  7  6  3  5  1  0  2 

CDAB  50317462 

BADC 

5  3  1  7  4  6  2  0 

DCBA  17250643 

CDAB 

which  is  an  0(12,2)  set. 
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Discussion.  The  necessary  requirements  for  the  construction  of  an  0(n,t)  set, 

n  =  nl  +  n2’  1  <  n2’  by  t*ie  method  of  sum  composition  are:  The  existence  of  an 

0(n^,t)  set,  n^>tn2,  with  at  least  tn^  common  parallel  transversals,  and  an 

0(n  ,t)  set.  These  conditions  are  obviously  satisfied  whenever  n,  and  n 
c  12 

are  prime  powers. 


While  for  some  values  of  n  there  exists  only  a  unique  decomposition  ful¬ 
filling  the  above  requirements,  for  infinitely  many  other  values  of  n  there  are 
abundant  such  decompositions. 

It  seems  that  if  there  exists  an  0(n.,2)  set  and  if  n  =  n.  +  n  ,  n.  >  2n_ 
then  one  can  construct  an  0(n,2)  set  by  the  method  of  sum  composition  if  n^ 
is  a  prime  power.  To  support  this  observation  and  shed  some  more  light  on  the 
method  of  sum  composition  we  present  in  subsequent  pages  some  highlights  of  the 


results  which  we  hope  to  complete  and  submit  for  publication  shortly. 


In  the  following  for  each  given  decomposition  of  n  we  exhibit  an  0(n,2) 


set  which  has  been  derived  by  the  method  of  sum  composition.  We  shall  represent 


the  pairs  in  a  form  that  the  curious  reader  can  easily  reconstruct  the  original  sets. 
Hereafter  the  notation  _l  means  that  is  orthogonal  to  . 

1)  12  =  9  +  3 


ABC 

4  5 

6 

7 

8 

9 

1 

2 

3 

12  3 

4  5 

6  A  B 

C 

8 

9 

7 

B  C  A 

1  2 

3 

4 

5 

6 

9 

7 

8 

9  7  8 

3  1 

2  B  C 

A 

6 

4 

5 

CAB 

7  8 

9 

1 

2 

3 

5 

6 

4 

5  6  4 

6  9 

7  C  A 

B 

1 

2 

3 

2  3  1 

5  6 

4 

A 

B 

c 

8 

9 

7 

6  4  5 

A  B 

C  3  1 

2 

7 

8 

9 

8  9  7 

2  3 

1 

B 

C 

A 

4 

5 

6 

2  3  1 

B  C 

A  8  9 

7 

5 

6 

4 

5  6  A 

8  9 

7 

C 

A 

B 

3 

1 

2 

x 

7  8  9 

C  A 

B  4  5 

6 

3 

1 

2 

3  12 

A  B 

C 

9 

7 

8 

6 

4 

5 

ABC 

2  3 

15  6 

4 

9 

7 

8 

9  7  8 

B  C 

A 

6 

4 

5 

2 

3 

1 

B  C  A 

7  8 

9  12 

3 

4 

5 

6 

6  4  5 

C  A 

B 

3 

1 

2 

7 

8 

9 

CAB 

6  4 

5  9  7 

8 

2 

3 

1 

1  5  9  6  7 

2 

8 

3  4 

A  B 

C 

4  9  2  5  7 

3  6  8  1 

A  B 

C 

7  2  6 

3  4  8 

5 

9 

1 

B 

C  A 

3  5  7 

1  6  8  2  4 

9 

CAB 

4  8  3 

9  1 

5 

2 

6  7 

CAB 

8  16 

9  2  4  7  3 

5 

B 

C  A 
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2)  14  =  11  +  3,  the  only  decomposition  which  fulfills  the  necessary  require¬ 

ments. 


A 

B 

C 

3 

4 

5 

6 

7 

8 

9 

10 

0 

1 

2 

B 

C 

9 

10 

0 

1 

2 

3 

4 

5 

A 

6 

7 

8 

C 

4 

5 

6 

7 

8 

9 

10 

0 

A 

B 

1 

2 

3 

10 

0 

1 

2 

3 

4 

5 

6 

A 

B 

C 

7 

8 

9 

6 

7 

8 

9 

10 

0 

1 

A 

B 

C 

5 

2 

3 

4 

2 

3 

4 

5 

6 

7 

A 

B 

C 

0 

1 

8 

9 

10 

9 

10 

0 

1 

2 

A 

B 

C 

6 

7 

8 

3 

4 

5 

5 

6 

7 

8 

A 

B 

C 

1 

2 

3 

4 

9 

10 

0 

1 

2 

3 

A 

B 

C 

7 

8 

9 

10 

0 

4 

5 

6 

8 

9 

A 

B 

C 

2 

3 

4 

5 

6 

7 

10 

0 

1 

4 

A 

B 

C 

8 

9 

10. 

0 

1 

2 

3 

5 

6 

_7 

0 

5 

10 

4 

9 

3 

8 

2 

7 

1 

6 

A 

B 

C 

7 

1 

6 

0 

5 

10 

4 

9 

3 

8 

2 

B 

C 

A 

3 

8 

2 

7 

1 

6 

0 

5 

10 

4 

9 

C 

A 

B 

0 

1 

2 

3 

4 

5 

6 

7 

A 

B 

C 

9 

10 

8 

8 

9 

10 

0 

1 

2 

3 

A 

B 

C 

7 

5 

6 

4 

5 

6 

7 

8 

1 

10 

A 

B 

C 

3 

4 

1 

2 

0 

2 

3 

4 

5 

6 

A 

B 

C 

10 

0 

1 

8 

9 

7 

10 

0 

1 

2 

A 

B 

C 

6 

7 

8 

9 

4 

5 

3 

7 

8 

9 

A 

B 

C 

2 

3 

4 

5 

6 

0 

1 

10 

4 

5 

A 

B 

C 

9 

10 

0 

1 

2 

3 

7 

8 

6 

1 

A 

B 

C 

5 

6 

7 

8 

9 

10 

0 

3 

4 

2 

A 

B 

C 

1 

2 

3 

4 

5 

6 

7 

8 

10 

0 

9 

B 

C 

8 

9 

10 

0 

1 

2 

3 

4 

A 

6 

7 

5 

C 

4 

5 

6 

7 

8 

9 

10 

0 

A 

B 

2 

3 

1 

6 

10 

3 

7 

0 

4 

8 

1 

5 

9 

2 

A 

B 

C 

3 

7 

0 

4 

8 

1 

5 

9 

2 

6 

10 

C 

A 

B 

9 

2 

6 

10 

3 

7 

0 

4 

8 

1 

5 

B 

C 

A 

3)  15  =  12  +  3,  15  =  11  +  4  are  the  only  decompositions  which  fulfill  the  neces¬ 

sary  requirements.  However,  we  consider  here  the  latter  decomposition  since 
we  can  utilize  the  properties  of  Galois  field  GF(ll). 
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1?  =  13  +  4  and  17  =  12  +  5  are  the  only  decompositions  which  fulfill  the  necessary  requirements. 
The  following  pair  is  derived  through  the  first  decomposition. 
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We  do  not  know  if  an  0(14,2)  set  with  8  common  parallel  transversals  exists  or  if  an  0(15,2) 
set  with  6  common  transversals  corresponding  to  the  0(15,3)  set  in  section  V  can  be  combined 
with  an  0(3,2)  set  to  form  an  0(18,2)  set.  Therefore,  the  following  0(18,2)  set  is  constructed 
from  the  decomposition  18  =  13  +  5  . 
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6)  We  do  not  know  if  an  0(15,2)  set  with  14  common  parallel  transversals 
corresponding  to  the  0(15,3)  set  in  section  V  can  be  combined  with  an 
0(7,2)  set  to  form  an  0(22,2)  set  or  if  there  exists  an  0(18,2)  set 
with  8  common  parallel  transversals.  Therefore,  the  following  0(22,2) 
sets  are  derived  from  the  decompositions  22  =  19  +  3  and  22  =  17  +  5  . 
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XIV.  Computer  Construction  of  0(10,  t)  Sets 


In  about  fifteen  years  the  effectiveness  of  computers  in  searching  for 
orthogonal  sets  of  latin  squares  of  order  ten  has  increased  strikingly.  Still 
the  problem  is  so  large  that  there  seems  to  be  little  reason  for  optimism  that  the 
order  ten  problem  can  be  completed  by  computers.  More  precisely,  if  (as  most 
conversant  with  the  problem  consider  quite  plausible)  no  0(10,3)  set  of 
orthogonal  latin  squares  of  order  ten  exists,  then  the  number  of  cases  to  con¬ 
sider  seems  too  large  for  an  exhaustive  proof  by  computer  to  be  achievable.  The 
number  of  latin  squares  of  order  ten  is  astronomical. 

About  19  53  Paige  and  Tompkins  [I960]  programmed  SWAC  to  search  for 
squares  orthogonal  to  a  fixed  latin  square  of  order  ten.  A  few  hours  of  running 
produced  no  orthogonal  square,  and  was  regarded  as  a  bit  of  experimental  evi¬ 
dence  for  the  truth  of  Euler's  conjecture.  Calculations  based  on  the  progress 
made  in  the  search  led  to  the  extrapolation  that  over  fifty  million  years  of  computer 
time  would  be  required  to  search  for  all  squares  orthogonal  to  a  latin  square  of 
order  ten  put  into  SWAC  intially.  (At  about  the  same  time  a  similar  program  was 
written  and  similar  results  obtained  with  MANIAC  at  Los  Alamos;  this  attempt 
has  not  been  reported  in  print. ) 

In  19  59,  after  Euler's  conjecture  had  been  disproved  for  all  orders 
41  +  2  >  6,  Parker  programmed  UNIVAC  1206  to  search  for  squares  orthogonal 
to  a  latin  square  of  order  ten.  The  running  time  was  sharply  less  than  for  SWAC 
or  MANIAC,  about  thirty  minutes  for  the  majority  of  latin  squares.  This  was 
accomplished  by  generating  and  storing  all  transversals  of  the  input  latin  square, 
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then  searching  for  all  ways  to  form  latin  squares  from  the  list  of  transversals. 

(A  transversal,  or  directrix,  is  a  set  of  cells  of  a  latin  square,  one  in  each  row, 
one  in  each  column,  and  one  containing  each  digit.)  The  striking  gain  in  speed 
over  the  earlier  efforts  occurred  largely  because  the  number  of  transversals  of 
a  typical  latin  square  of  order  ten  is  roughly  850,  much  less  than  10!;  and, 
of  course,  the  search  was  several  levels  deep.  (SWAC  and  MANIAC  were  pro¬ 
grammed  to  build  up  starts  of  latin  squares  toward  orthogonal  mates  by  filling  in 
cells  to  form  rows.) 

There  were  two  main  outcomes  from  considerable  running  of  Parker's  1206 
program:  1)  0(10,3)  sets  of  latin  squares  are  not  numerous;  more  precisely, 
only  a  small  fraction,  if  any,  order  ten  squares  could  possibly  extend  to  0(10,3) 
sets.  Some  400  latin  squares  were  run.  Some  were  random,  some  were  com¬ 
puter  output  fed  back  as  input  and  hence  known  to  have  an  orthogonal  mate,  and 
some  were  considered  interesting  candidates  for  intuitive  reasons  by  Parker  and 
others.  Not  once  did  an  exhaustive  search  for  orthogonal  mates  of  an  input  latin 

square  discover  a  pair  orthogonal  to  one  another.  Mild  evidence  may  be  claimed 
supporting  the  opinion  that  no  0(10,3)  set  exists.  2)  Of  a  computer-generated 

sample  of  100  random  latin  squares  of  order  ten  (program  by  R.  T.  Ostrowski), 

62  have  orthogonal  mates.  Thus,  unlike  0(10,3)  sets,  0(10,2)  sets  of  squares 
are  quite  common.  Euler's  intuition  for  order  ten  was  not  only  wrong,  but  in  this 
sense  wrong  by  a  large  margin.  It  was  this  finding  which  tempted  Parker  for  a 
time  to  believe  that  repeated  runs  of  the  program  should  have  a  good  chance  of 
producing  at  an  0(10,3)  set,  but  many  failures  dimmed  this  optimism. 
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In  1967  John  W.  Brown  programmed  IBM  7094  to  decide  whether  an  input 
latin  square  of  order  ten  can  be  extended  to  an  0(10,3)  set.  The  running  time 
was  one  half  minute.  Almost  needless  to  say,  transversals  again  were  generated. 
Searching  for  patterns  of  transversals  toward  extension  to  an  0(10,3)  set  pro¬ 
duced  a  speed  gain  over  the  previous  program  for  orthogonal  pairs.  Brown 
endeavored  to  get  every  drop  of  speed  from  the  machine.  As  before,  hundreds 
of  input  order-ten  latin  squares  produced  no  0(10,3)  set. 
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XV.  On  the  Equivalence  of  0(n,  t)  Sets  With  Other  Combinatorial  Systems 


XV.  0.  Summary 

In  this  section  we  have  densely  summarized  some  of  the  results  obtained 
by  author  and  at  least  fourteen  others  in  order  to  demonstrate  the  importance  of 
the  theory  of  mutually  orthogonal  latin  squares.  We  have  shown  that  fourteen 
well-known  and  important  combinatorial  systems  with  certain  parameters  are 
actually  equivalent  to  a  set  of  mutually  orthogonal  latin  squares.  A  schematic 
representation  of  these  equivalences  has  been  demonstrated  in  four  wheels  which 
we  have  called  "Fundamental  Wheels  of  Combinatorial  Mathematics". 

XV.  1.  Introduction 

The  theory  of  mutually  orthogonal  latin  squares  owes  its  importance  to  the 
fact  that  many  well-known  combinatorial  systems  are  actually  equivalent  to  a 
set  of  mutually  orthogonal  latin  squares;  viz.,  finite  projective  plane,  finite 
Euclidean  plane,  net,  BIB,  PBIB,  orthogonal  arrays ,  a  set  of  mutually  orthogonal 
matrices,  error  correcting  codes,  strongly  regular  graphs ,  complete  graphs,  a 
balanced  set  of  i-restrictional  lattice  designs,  difference  sets,  Hadamard 
matrices,  and  an  arrangement  of  non  attacking  rooks  on  hyperdimensional  chess 
board.  These  combinatorial  systems  are  unquestionably  potent  and  effective  in 
all  branches  of  combinatorial  mathematics,  and  in  particular,  in  the  construction 
of  experimental  designs.  Therefore,  a  statement  that  the  theory  of  mutually 

i 

t 

orthogonal  latin  squares  is  perhaps  the  most  important  theory  in  the  field  of 
experiment  designs  is  not  in  the  least  exaggerated  as  far  as  this  author  is  con¬ 
cerned. 
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Our  purpose  in  this  section  is  to  demonstrate  the  relation  of  a  set  of 
mutually  orthogonal  latin  squares  with  the  above  mentioned  combinatorial  systems. 
We  shall  present  the  essence  of  the  known  results  available  only  in  scattered 

4 

literature  in  one  theorem  which  we  consider  to  be  a  "fundamental  theorem  of 
combinatorial  mathematics".  For  the  definitions  of  these  combinatorial  systems 
and  the  proof  of  the  forthcoming  theorem  see  the  list  of  references  given  at  the 
end  of  this  paper. 

XV.  2.  Notation 

For  the  sake  of  conciseness  we  introduce  the  following  notations: 

0)  0(n,t)  denotes  a  set  of  t  mutually  orthogonal  latin  squares  of  order 

n  . 

,  1)  MOM(n,t)  denotes  a  set  of  t  mutually  orthogonal  n  X  n  matrices. 

2 

2)  OA(n,t)  denotes  a  set  of  orthogonal  arrays  of  size  n  ,  depth  t,  n 
levels,  and  strength  2  . 

3)  Net(n,t)  denotes  a  net  of  order  n  and  degree  t  . 

4)  Code(n,r,t;m)  denotes  a  set  of  n  code  words  each  of  length  r  such 
that  any  two  code  words  are  at  least  at  Hamming  distance  >  t  on  an 
m-set  2  with  m  distinct  elements.  We  remind  the  reader  that  such 
a  code  is  also  called  (t-l)-error  detecting  code  or  (t-l)/2-error  cor¬ 
recting  code  because  such  a  code  is  capable  of  detecting  up  to  t-1 
errors  and  correct  up  to  (t— 1 )/ 2  errors  in  each  transmitted  code  word. 

5)  PBIB(b, v,r,k denotes  a  partially  balanced  incomplete  block 

.  design  with  b  blocks  each  of  size  k,  v  treatments  with  r  replication 

of  each,  and  association  indices  and  . 
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6)  SR-Graph  (A)  denotes  the  strongly  regular  graph  with  incidence 
matrix  A  . 

7)  Non  #(n,t)  denotes  an  arrangement  of  n  mutually  non  attacking 
rooks  on  the  t-dimensional  n  X  n  chess  board. 

8)  PG(2,s)  denotes  a  finite  projective  plane  of  order  s  (not  necessarily 
Desarguesian). 

9)  £(2,s)  denotes  a  finite  Euclidean  plane  of  order  s  . 

10)  BIB(b,  v,r,k,\ )  denotes  a  balanced  incomplete  block  design  with  b 
blocks  each  of  size  k,  v  treatments  with  r  replications  of  each,  and 
association  index  \  . 

11)  K-Graph  (A)  denotes  the  complete  graph  with  incidence  matrix  A  . 

12)  DIF(v,k,\)  denotes  a  difference  set  with  parameters  v,  k,  and  X.  . 

13)  BLRL(s)  denotes  a  balanced  set  of  i-restrictional  lattice  design  for 
s  treatments.  Note  that  a  1-restrictional  balanced  lattice  design  is 
simply  a  BIB  design. 

14)  HAD(n)  denotes  a  symmetric  normalized  Hadamard  matrix  of  order  n  . 

Hereafter  we  also  adopt  the  following  two  notations: 

i)  A  <=>  B  means  A  implies  B  and  B  implies  A  . 

ii)  A  =>  B  means  A  implies  B  .  Whether  or  not  B  implies  A  is 

unstated. 
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XV.  3.  The  Result 


Theorem 

(a)  For  any  pair  of  positive  integers  n  and  t  we  have: 

1)  0(n,t)  <==>  MOM(n,t+2) 

2)  0(n,t)  <=s>  OA(n,t+2) 

3)  0(n,t)  <=>  Net(n,t+2) 

2 

4)  0(n,t)  <*=>  Code(n  ,t+2,  t+l;n) 

5)  0(n,t)  PBIB(n2,n(t+2),  n,t+2,0,l) 

6)  0(n,t)  <=£>  SR-Graph  (A)  where  A  is  the  incidence  matrix  associated 

with  PBIB  in  5). 

7)  0(n,t)  <==>  Non  #(n2,nt+2)  . 

(b)  H  t  =  n-1  then  also: 

8)  0(n,n-l)  »  PG(2,n) 

9)  0(n,n-l)  <=>  £(2,n) 

2  2 

10)  0(n,n-l)  <=>  BIB(n  +n+l,  n  +n+l,  n+1,  n+1,1) 

2  2 

11)  0(n,n-l)  <=£>  Code(n  +n+l,  n  +n+l,  2n;2) 

12)  0(n,n-l)  K-Graph  (A)  where  A  is  the  incidence 

matrix  associated  with  BIB  in  10) 

13)  0(n,n-l)  <=£>  DIF(n2+n+l,  n+1,1)  . 

(c)  If  n  =  pm  where  p  is  a  prime  and  m  is  a  positive  integer  then  also  the 
following; 

14)  0(pm,pm-l)  <=>  BLRL(pm)  . 
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(d)  If  n  =  2r  and  t  =  r- 2,  r  >  3  then  the  following  are  also  true: 


15) 

0(2r,r-2) 

=i>HAD(4r2) 

16) 

0(2r,  r-2) 

==>  BIB(4r2  -  1,  4r2-l,  2r2-l,  2r2-l,  r2-l) 

17) 

18) 

0(2r,  r-2) 

0(2r,r-2) 

=>  Code(4r2-l,  4r2-l,  2r2;2) 

=>  Code(8r2,  4r2,  2r2;2) 

19) 

0(2r ,  r-2) 

=£>  DIF(4r2-l,  2r2-l,  r2-l)  . 

A  complete  schematic  representation  of  this  theorem  can  be  demonstrated 
in  four  wheels  which  will  be  called  "fundamental  wheels  of  combinatorial  mathe 
matics".  For  the  sake  of  compactness  we  shall  omit  the  associated  parameters 
with  each  system  in  these  wheels  except  for  0(n,t)  .  By  knowing  the  values 
of  n  and  t  in  the  given  0(n,t)  sets, then  the  reader  can  easily  find  the 
associated  parameters  with  other  systems  in  the  wheels  from  the  proper  part  of 
the  above  theorem. 
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Wheel  1.  For  any  positive  integer  n  and  t 
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Wheel  2.  For  any  positive  integer  n 
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Wheel  3. 


For  any  prime  p 


and  positive  integer 


m  . 
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